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ABSTRACT 


Sea  surface  motion  has  been  observed  for  two  decades  by  viewing  successive  satellite 
thermal  infrared  Images  for  movement  of  thermal  gradients.  A  major  effort  has  been  made 
over  the  last  eight  years  to  automate  the  task  of  measuring  surface  current  flow  using 
satellite  images,  a  task  made  difficult  by  the  presence  of  land,  clouds  and  other  unwanted 
features. 

This  paper  preseats  a  method  to  quantitatively  measure  ocean  surface  movement. 
Using  10.5  fim  band  imagery,  sequential  images  are  processed,  correcting  for  the  presence 
of  clouds  and  land,  to  determine  flow  direction,  magnitude,  and  velocity.  An  ordered 
statistical  edge  detection  algorithm  is  used  to  detect  gradient  boundaries,  discriminate 
between  the  water  surface,  land  and  clouds,  and  select  thermal  features.  A  constrained 
correlation  based  pattern  recognition  scheme  is  then  used  to  match  the  feature  in  the  next 
image.  Surface  flow  is  calculated  for  each  selected  point  with  results  presented  in  tabular 
and  image  form.  A  study  area  in  the  region  of  the  Delaware  Bay,  on  the  east  coast  of 
the  United  States,  has  been  analyzed  using  this  technique,  with  the  results  comparing 
favorably  with  in  situ  anchored  buoy  measurements. 


Chapter  1 


INTRODUCTION 


1.1  Purpose 

Ima^ng  satellites  provide  viewing  of  the  earth’s  surface  in  temporal  and  spatial 
scales  that  cannot  be  matched  by  surface  bound  collection  methods.  Modeling  of  the 
world’s  oceans  mixed  layer  depths,  a  vital  ocean  region  for  primary  production  of  food 
and  oxygen,  is  very  dependent  on  surface  current  information.  At  present,  most  of  this 
information  is  available  only  from  ships,  research  vessels  and  anchored  or  free  floating 
buoys.  These  platforms  have  small  observation  areas,  providing  only  the  most  limited 
coverage  of  the  ocean  surface. 

This  work  was  initiated  to  provide  the  Center  for  Remote  Sensing,  University  of 
Delaware,  with  a  method  to  objectively  process  satellite  thermal  images  to  estimate  sea 
surface  velocity. 

1.2  Use  of  satellites  to  view  the  earth  surface 

Images  collected  and  sent  by  satellites  offer  unprecedented  viewing  of  the  planet  we 
live  on.  Satellites  in  orbit  today  for  civil  applications  are  viewing  the  earth  with  visible, 
near  infrared,  thermal  infrared  and  radar  sensors.  Figure  1.1  illustrates  the  key  concepts 
of  satellite  imaging  of  the  earth’s  surface. 

Geostationary  (GOES)  satellites  are  placed  in  a  very  high  orbit  above  the  earth 
(approximately  36,000  km;  to  have  their  orbit  match  the  earth’s  rotation  rate.  This  allows 
viewing  of  the  same  area  at  all  times,  but  limits  resolution  to  very  large  features. 

Low  orbiting  satellites  (500-2,000  km)  travel  around  the  earth  at  velocities  relating 
to  their  altitude.  Most  mapping  satellites,  including  the  ones  that  provided  the  images 
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2 

GOES  (Geostationaiy) 
I  '"'"I  Satellite 


Figure  1.1:  Low  orbiting  and  geostationary  satellite  comparison. 


for  this  study,  are  in  low  polar  orbits  where  the  earths  rotation  moves  perpendicular  to 
the  satellite  path,  allowing  coverage  of  a  portion  or  the  entire  surface  over  time. 

1.3  Satellite  imagery  used  in  this  study 

This  paper  uses  images  from  National  Oceanic  and  Atmospheric  Administration 
(NOAA)  low  orbiting  satellites,  NOAA  -  8,  9  and  10.  The  images  used  in  this  study  are 
from  the  10.5  fim  band.  Advanced  Very  High  Resolution  Radiometer  (AVHRR)  sensor. 
They  are  in  a  512x512  pixel  format,  and  are  geometric  and  grey  scale  value  for  temperature 
corrected,  with  a  nominal  pixel  dimension  of  1.1  kilometers  at  the  nadir.  From  a  series  of 
images  covering  the  period  1  January  to  30  June  1989,  sets  acceptably  free  from  clouds 
were  selected  for  the  following  periods: 


•  10-11  March  1989. 
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•  28-29  May  1989. 

•  10-12  June  1989. 

•  29-30  June  1989. 

The  images  are  stored  in  8  bit,  0-255  value,  grey  scale  form.  Ocean  thermal  features 
are  visible  as  gradients  (Figure  1.3).  Temperature  can  be  determined  from  the  image  using 
a  linear  function  with  a  slope  of  0.125  and  a  y  intercept  of  zero  in  the  form: 

T  =  V'p.>e/X  0.125.  (1.1) 

This  gives  the  image  a  scaling  temperature  range  from  0  -  32°  Centigrade  (C).  Features 
with  surface  temperatures  above  or  below  this  range  are  treated  as  sensor  saturating 
signals  and  are  assigned  255  or  0  respectively  (Figure  1.2). 

1.4  Principles  of  ocean  surface  flow  estimation 

The  work  of  Vastano  and  Reid  [1],  categorizing  surface  flow  in  the  Oyashio  Frontal 
Zone  in  the  northwest  Pacific  Ocean  was  the  starting  point  for  this  project.  Their  ob¬ 
servations  using  AVHRR  (ll/im  band)  images  showed  an  11°  centigrade  (C)  temperature 
variation  in  Sea  Surface  Temperature  (SST).  They  observed  defined  regions  of  tempera¬ 
ture  which  they  associated  with  frontal  regions  of  differing  streamfiow.  These  temperature 
gradients  are  produced  by  ‘^small  scale  waves  or  turbulent  processes  ...  of  (what  should 
be)  otherwise  smooth  surface  isotherms.” 

Rossby  like  internal  gravity  waves,  the  effects  of  heat  exchange  with  the  atmosphere, 
and  vertical  motion  on  the  sea  surface  were  assumed  to  produce  much  smaller  changes  in 
SST  than  the  turbulance  caused  by  currents  which  bring  core  temperatures  to  the  surface. 

1.4.1  Surface  flow  model 

Rossby  like  waves  have  a  length  and  time  scale  described  by  [2]: 

r  _  ^^0 

~h' 
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Figure  1.2:  NOAA  thermal  image  temperature  scale.  One  pixel  value  equals  a  change 
of  0.125  degrees  centigrade. 

and 

_  2k  R  tanO 

T  = - 

NoD  ’ 

where  L  is  the  Rossby  scale  of  horizontal  motion,  No  is  the  characteristic  buoyancy,  D  is 
the  vertical  scale  characteristic,  fo  is  the  coriolis  parameter,  R  is  the  radius  of  the  earth, 
and  0  is  the  latitude. 

Using  values  from  Stewart  [2],  for  40®  North  (the  region  of  interest  for  the  images 
analyzed  in  Vastano’s  study)  yields: 


X  =  74  km. 


Figure  1.3:  Typical  satellite  thermal  image  of  coastal  Delaware  and  the  continental  shelf. 

0132  hours  local,  30  June  1989  (false  color  added  for  feature  enhancement). 

r  =  83  days, 

Flownosiby  =  0.011  m/s. 

For  study  of  the  Delaware  -  New  Jersey  coastal  area  (30°  North),  with  measured 
flows  of  0.1  to  1.2  meters  per  second,  Vcistano’s  assumption  to  ignore  Rossby  waves  still 
applies,  since  Flowngaaby  «  0.013  m/s. 

In  Vastano  and  Reid’s  paper,  local  changes  in  sea  surface  temperature  are  described 
with  the  following  model: 


where  C  represents  atmosphere  heat  exchange,  diffusion  and  vertical  motion  perturbation 
components,  V  is  the  perturbation  due  to  flow  and  R  is  the  rate  of  temperature  change. 
C  and  V  each  have  two  components: 


C  =  Co  +  Cl , 

(1.3) 

V  =  Vo+Vi, 

(1.4) 

where  Cq  and  Vb  are  mean  components  of  change  and  C\  and  Vi  represent  the  perturbation 
effects. 

Expanding  equation  1.2  with  1.3  and  1.4  yields: 


aC'i  =  7?-Vb-VCi  =  5i  +  S2,  (1.5) 

where  S\  is  the  self  advective,  turbulent  diffusion  term.  Expanding  Si  to  represent  effects 
of  wind,  convection  and  internal  waves: 


Si  =  iZ  -  Vo  •  VCo  -  V, .  VCi  -  ^  . 

at 


S2  is  the  temperature  change  due  to  flow: 


(1.6) 


S2  =  -Vi.VCo,  (1.7) 

Using  the  observation  that  the  mean  terms  of  C  and  V  change  slowly,  relative  to 
the  rapid  changes  of  temperature  over  distance  and  time  on  the  flow  boundaries,  equation 
1.6  is  disregarded.  This  leaves  equation  1.7  as  the  representing  source  term  to  describe 
changes  in  Sea  Surface  Temperature. 

Away  from  flow  boundaries  this  change  can  be  described  by: 
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where  the  velocity  Vq  represents  the  velocity  of  small  discrete  volumes  of  water.  This 
velocity  is  assumed  to  be  equivalent  to  the  time  rate  of  change  of  the  displacement  of  the 
parcel  between  sampling  periods  (images). 

Assuming  the  flow  is  horizontally  nondivergent  and  using  an  approximation  of  a 
uniform  Coriolis  parameter,  Vastano  concluded  that  the  surface  velocity  can  be  approxi¬ 
mated  by  a  single  scalar  streamfunction,  This  function  is  represented  in  trigonometric 
terms  by: 


_  .  .  nttx  .  miry 

^  =  2^^  ■^n.m  sin  sin 

n  m  * 


Ly  ’ 


(1.9) 


where  Lx  and  Ly  are  the  dimensions  of  the  region  surrounding  the  samples. 

The  vertical  («)  and  horizontal  (t>)  velocity  components  are  derived  from: 


u  = 


dy 


.  nirx  mir 

sin  — —  cos  , 

Lx  Ly 


(1.10) 


nir 


An, 


nnx  mir 
,nCOS-=— cos  — 


(1.11) 


1.4.2  Additional  considerations  for  the  Vastano  and  Reid  model 

Wahl  and  Simpson  [3],  have  shown  that  the  assumptions  of  the  Vastano  and  Reid 
model  concerning  horizontal  dilTusion  and  changes  in  Sea  State  Temperature  (SST)  are  in 
error.  At  wind  surface  wind  speeds  less  than  10  meters /second  (m/s),  horizontal  as  well 
as  vertical  diffusion,  and  air  -  sea  heat  exchange  are  the  dominant  process  in  SST  change. 
Surface  winds  greater  than  20  m/s  usually  cause  intense  vertical  mixing,  destroying  the 
thermal  feature.  Work  by  Price  [4]  has  shown  that  diurnal  heating  also  causes  significant 
distortion  of  thermal  features  by  wind  stress  induced  diurnal  jets  which  transport  surface 
thermal  features  normal  to  wind  direction  at  velocities  up  to  0.1  m/s  for  moderate  winds. 

The  conservation  of  heat  equation  now  becomes: 
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—  =  iir»v5,r  +  jir,^  +  5;Oi,  (1.12) 

t 

where  ^  is  the  material  derivative  of  temperature  over  time,  Kh  is  the  horizontal  diffu- 
sivity,  Qi  is  the  heating  source  (or  sink)  term,  T  is  the  temperature,  an  t  is  time.  The 
derivative  of  temperature  is  defined  by: 


dT  dT 

■5r  =  '5r+“ 


VT, 


(1.13) 


where  tt  =  (u,  v,  w)  is  the  fluid  velocity.  For  a  two  -  dimensional,  non  -  compressible  flow, 
assuming  a  well  mixed  layer  (U  <  |f)  and  |^),  ^  becomes: 


—  ^U^VT  +  KHVjjT  +  Y^Qi,  (1.14) 

t 

where  the  first  term,  U  •  VT,  is  the  advection  term,  the  second  term,  KjfV^T,  is  the 
horizontal  diflTusion  term  and  the  last  term,  Qi  represents  the  air  -  sea  interaction. 


1.5  The  validity  of  surface  feature  duration  assumptions 

The  fundamental  assumption  which  allows  estimation  of  surface  flow  from  satellite 
imagery  is  that  the  thermal  feature,  representing  the  temperature  of  a  volume  of  water, 
remains  relatively  intact  during  the  sampling  interval.  Each  pixel  represents  1.21  square 
kilometers  (km^)  or  greater  of  sea  surface.  The  latent  heat  flux  calculated  using  satellite 
data  represents  the  majority  of  the  temperature  of  a  layer  of  homogeneous  temperature 
water,  called  the  surface  mixed  layer.  This  mixed  layer  may  extend  to  depths  of  a  few 
hundred  meters  in  winter,  to  less  than  40  meters  in  summer. 

The  critical  assumption  which  allows  satellite  determination  of  surface  flow  is  that 
thermal  features  sampled  over  a  large  area  (  a  10  x  10  pixel  sample  is  an  area  greater  than 
127  km?  ),  retain  their  three  dimensional  identity  over  periods  of  up  to  24  hours.  This 
assumption  corresponds  to  one  consecutive  pass  by  a  single  low  orbit  satellite  over  the 
same  location.  This  is  true  as  long  as  surface  flow  is  dominated  by  by  geostropic  currents, 
resulting  in  only  advective  motion  in  the  horizontal  plane,  with  little  or  no  horizontal 
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diffusion.  The  presence  of  wind  driven  diurnal  currents,  giving  Kjj  velocities  of  5  -  10 
cm/s  casts  some  question  as  to  the  validity  of  the  assumption  that  horizontal  diffusion  is 
minor  term.  Answering  this  question  is  beyond  the  scope  of  this  work.  The  assumption 
A'^V^T  d  U  •  VT  is  considered  valid  for  periods  less  than  24  hours. 


Chapter  2 


THERMAL  FEATURE  SELECTION  AND  TRACKING 


2.1  Feature  selection 

Using  the  assumption  that  a  large  volume  of  water  will  hold  a  distinctive  thermal 
signal  shape  over  time,  it  is  possible  to  employ  pattern  recognition  techniques  to  measure 
displacement  of  features  defined  by  temperature  alone.  Once  features  are  selected  from  an 
initial,  or  pattern  image,  motion  can  be  estimated  by  locating  them  and  measuring  their 
displacement  in  subsequent  images.  In  this  work,  this  accomplished  in  three  steps:  Pattern 
features  are  first  selected  from  an  initial  image  using  Ordered  Statistical  Edge  Detection  or, 
for  comparison  to  other  models,  grid  sampling  methods.  Use  of  edge  detection  to  extract 
features  insures  that  only  areas  of  gradient  change  are  mapped,  eliminating  samples  of 
featureless  areas  that  grid  methods  produce.  Second,  feature  recognition,  classification 
and  selection  are  performed  using  a  comparison  image,  searching  areas  having  maximum 
likelihood  of  containing  the  matching  pattern.  As  the  last  step,  displacement,  direction 
and  velocity  calculations  are  made  based  on  movement  of  the  matched  features  relative 
to  the  pattern  features  using  a  common  coordinate  system. 

Edge  detection  (Figure  2.1)  was  employed  to  select  features  that  are  involved  in 
advective  flow  across  the  sea  surface,  as  opposed  to  more  standard  grid  methods  used  by 
Kelly  [5]  and  Simpson  [6],  to  locate  areas  of  maximum  flow.  This  approach  reduces  the 
number  of  points  selected,  with  a  corresponding  reduction  in  the  number  of  calculations 
required  to  process  an  image.  Figure  2.2  shows  a  three  -  dimensioned  representation  of  the 
thermal  signal  for  the  Delaware  Bay  mouth  and  coastal  area,  one  of  the  areas  of  interest  to 
this  study.  Feature  selection  maps  an  area  of  this  signal,  discriminating  against  selection 
of  land,  clouds  and  other  unwanted  features. 
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Figure  2.1:  Thermal  features  detected  by  the  SEDA  method,  superimposed  on  Figure 
1.3,  using  a  5  X  5  edge  detection  tile  {Threshold  =  6,  Range  =  12). 

2,1.1  Ordered  statistical  edge  detection 

The  Ordered  Statistical  Edge  Detection  Algorithm  (SEDA)  which  is  used  for  feature 
selection  is  one  adapted  from  the  work  of  Hardie  and  Gnacek  [7],  based  on  the  paper  by 
Pitas  and  Venetsanopoulos  [8].  This  type  of  edge  detector  has  proven  to  be  very  effective 
in  noisy  signal  environments.  The  SEDA  method  is  not  orientation  dependent.  Even 
though  it  scans  from  left  to  right,  it  accurately  detects  ascending  and  descending  gradients, 
marking  curved  and  linear  features.  The  method  is  very  effective  in  locating  features  in 
noisy  signal  environments.  It  works  by  incrementally  moving  an  n  x  n  tile  across  an  image 
looking  for  thermal  gradients  shown  by  contouring  in  Figure  2.3. 

Using  the  SEDA  method,  each  pixel  is  defined  by  a  vector  in  a  set  of  n  vectors 


12 


Cape  May,  New  Jersey 


Thermal  map:  arbitrary  scale,  Azimuth:  3 IS 

Figure  2.2:  Three  dimensional  thermal  signal  map,  1258  hours  local,  30  June  1989, 
Delaware  Bay  and  coast.  The  high  values  correspond  to  land  temperatures 
for  Cape  May,  New  Jersey,  and  the  coast  near  Lewes,  Delaware.  The  very 
low  values  in  the  lower  right  corner  are  light  cumulus  clouds.  The  water 
temperature  signal  varies  considerably  from  the  Delaware  Bay  (upper  left, 
behind  Cape  May),  to  the  ocean  waters  in  the  foreground.  The  effects  of 
turbulent  mixing  can  be  easily  seen  in  the  bay  and  its  mouth  as  noticeable 
perturbations  in  surface  signal. 
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Cape  May,  New  Jersey 


10  20  30  40  50 

Thermal  Comour:  arbitrary  scale 


Figure  2.3:  Thermal  signal  contour  map,  Delaware  Bay  and  coastal  region,  1258  hours 
local  30  June  1989.  Cape  May  is  the  center  top.  Coastal  Delaware  is  on  the 
lower  left.  Light  clouds  can  be  observed  on  the  lower  right. 


defined  by  the  image  dimensions  such  that; 


Wn  —  ®1, 02>  •••>  ®n  • 

Each  vector  sample  is  a  set  of  n  vectors  defined  by  the  edge  detection  tile  dimensions: 

Q  =  [dl ,  02, ...,  On]  • 

The  distance  or  value  of  the  vector  is  the  pixel  value  (0  -  255): 


dn  =  llonll . 
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As  the  edge  detection  tile  moves  across  its  search  area,  it  generates  a  set  of  vectors: 


®2*  •••>  ®n  j 


with  values: 


^2»  •••>  • 

After  rank  ordering  from  highest  to  lowest,  the  set  notation  changes  to: 


0(1)  <  0(2)  <  —  <  0(„)  , 


with  vector  having  the  largest  value  of  d. 

The  SEDA  algorithm  selects  an  edge,  S,  when  the  range  exceeds  a  specified  value, 
called  a  Threshold  (T),  represented  in  equation  form  by: 


with  the  range  defined  as: 


«(•)  = 


I 


1 

0 


if  range  >  T 
otherwise 


(2.1) 


range  =  fl(n-i)  -  0(i+,) . 

The  value  t  is  the  number  of  ordered  statistics  chosen  toward  the  mean  to  evaluate  the 
sample  and  n  is  the  number  of  pixels.  The  range  value  is  the  minimum  change  in  intensity 
(pixel  value)  required  to  define  the  feature  as  an  edge.  Allowing  t  =  1,  for  example,  selects 
edges  based  on  n  —  2  samples,  ignoring  the  two  opposite  extreme  values.  The  feature  will 
be  recorded  if  the  0(„_q  —  0(i^q  value  is  greater  than  the  specified  T.  The  value  of  t  used 
determines  the  type  of  edge  selected.  As  i  approaches  n/2,  the  edge  selected  increases  in 
relative  sharpness,  defined  as  an  increasingly  steep  slope  in  the  linear  expression  of  the 
feature  pixel  values.  Since  the  range  criteria  must  be  met  with  a  difference  calculated 
from  a  decreasing  set  of  points,  the  slope  must  increase.  A  small  value  of  i  compared  to 
n  will  select  smoother,  less  sharp  edges.  A  value  of  f  =  1  was  used  in  this  study. 
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2.1.2  SEDA  feature  discrimination 

The  SEDA  algorithm  was  modified  during  this  project  to  allow  it  to  discriminate 
between  ocean  thermal  gradient  patterns  and  undesirable  features  such  as  cloud  and  land 
masses.  These  methods  are  not  one  hundred  percent  effective,  but  do  significantly  reduce 
feature  selection  error.  Figure  2.4  shows  pattern  tiles  mapped  for  a  study  of  the  Delaware 
coastal  region  selected  using  the  edge  detection  method.  The  method  selected  water 
features,  except  for  tile  00  wliich  was  chosen  in  an  area  of  light  clouds. 


Figure  2.4:  Pattern  tile  selection  in  the  Delaware  Bay  and  coastal  region,  1303  hours 
local.  10  June  19S9,  using  the  SED.A.  method.  One  falso  detection,  tile  00. 
north  of  the  bay  was  caused  by  light  clouds.  The  scale  in  the  lower  right 
is  the  distance  in  kilometers.  The  arrow  represents  displacement  for  a  flow 
rate  of  one  meter  per  second. 


The  modification  places  additional  restrictions  to  equation  2.1  to  take  adrantage 
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Table  2.1:  Typical  values  for  feature  selection  classification. 


Operator 

Value 

Use 

range 

6 

Minimum  Threshold:  change  must  be  greater  than  0.75®  C. 

rangemax 

12 

Maximum  Threshold:  change  must  not  be  greater  than  1.5®  C. 

maski 

131 

Minimum  Pixel  Value:  lowest  ranked  pixel  must  be  greater 
than  this  value  (summer  »  16.4®  C). 

mask2 

177 

Maximum  Pixel  Value:  highest  ranked  pixel  must  be  less  than 
this  value  (summer  »  22.13®  C). 

of  the  large  differences  in  pixel  values  between  the  ocean  and  other  features.  The  new 
selection  algorithm  is  described  by: 


1  if  T  <  range  <  rangcmax 
&  0(1+,)  >  maski 

&  a(n-«)  <  mask2 

0  otherwise 


(2.2) 


where  range^nax  specifies  a  maximum  change  that  can  occur  across  a  sample,  maski  is 
the  minimum  allowed  pixel  value,  and  mask2  is  the  maximum  allowed  pixel  value. 

Maski  and  mask2  are  determined  from  sampling  water  temperatures  in  the  area 
of  interest  and  applying  an  allowed  deviation  factor.  Range^ax  and  range  provide  a 
constraint  to  the  maximum  allowed  change  in  apparent  temperature.  Table  2.1  shows 
some  typical  values. 

2.1.3  Feature  selection  search  and  recording 

Once  features  are  selected  and  mapped  from  an  image,  the  features  are  referred  to 
as  Pattern  Tiles.  They  are  selected  by  searching  an  image  or  an  image  subset,  using  the 
modified  SEDA  algorithm  or  by  selecting  a  grid  of  points  to  be  sampled. 

1 

2. 1.3.1  SEDA  method 

The  SEDA  method  starts  by  evaluating  data  from  a  5  x  5  edge  detection  tile  (edge 
tile)  (Figure  2.5).  When  an  edge  tile  meets  the  modified  SEDA  criteria,  a  pattern  tile 
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maps  an  expanded  area  centered  on  the  edge  tile.  The  pixel  values  are  rank  ordered  and 
compared  with  the  mask  values.  If  the  high  and  low  extreme  mapped  values  are  within 
the  mask  values,  the  pattern  tile  is  stored,  the  counter  is  incremented  and  a  variable  offset 
is  applied  to  allow  spacing  between  the  selected  areas.  The  selection  continues  until  the 
search  area  specified  is  exhausted.  For  this  study,  an  11  x  11  pattern  tile  was  used. 


Figure  2.5:  SEDA  feature  selection  and  discrimination.  Selected  tiles  are  on  gradients 
in  regions  of  moderate  temperature  change,  free  from  clouds  and  land.  Tiles 
not  selected  are  rejected  if  the  temperature  in  the  mapped  area  exceeds 
classification  criteria. 


2. 1.3.2  Grid  search  method 

When  regular  spacing  of  points  is  desired,  a  grid  pattern  can  be  specified.  Entering 
the  coordinates  of  the  pattern  area,  the  program  will  sample  an  11  x  11  tile  at  each  point. 
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If  the  tile  is  acceptable,  it  is  stored,  the  counter  is  incremented  and  the  selection  continues 
until  the  area  is  exhausted. 


2.2  Feature  tracking 
2.2.1  Search  area  definition 

Once  the  initial  image  has  been  processed  and  pattern  tiles  have  been  selected,  the 
second  image  is  used  for  feature  tracking.  The  size  of  the  area  to  be  searched  for  each 
selected  feature  is  governed  essentially  by  the  Bayes  Decision  rule  [9]: 


P(6\C)  = 


PiC  I  0)Pi») 


where  P{6  \  ^  is  the  a  posteriori  probability  that  an  acceptable  match  will  be  found  given 
an  area  to  search  defined  by  the  vector  C*  The  likelihood  that  an  acceptable  match  can 
be  found  is  expressed  by: 


^(Cl  = 


fi(Cl 
^2(Cl  ’ 


where  Pi(C}  probability  that  the  feature  is  detectable,  and  lies  in  the  search  area. 
P2(^  is  the  probability  that  the  feature  can  be  recognized  and  selected  somewhere  in  the 
image.  Equations  2.3  and  2.4  can  be  expressed  as  a  maximum  a  posteriori  probability 
(MAP)  using  the  minus  •  log  likelihood  expression: 


^lnP{0  I  cl  =  ^InPiC  I  0)  +  §^lnP{0) ,  (2.5) 

where  the  first  term  on  the  right  represents  a  priori  information,  and  the  second  term 
is  a  posteriori  information.  Knowing  that  the  selected  features  are  altered  over  time  by 
diffusion  and  rotation,  and  the  displacement  velocity  varies  greatly  over  short  distances, 
equation  2.5  can  only  be  solved  through  nonparametric  density  estimation.  This  is  difficult, 
requiring  a  posteriori  information.  The  estimate  would  certainly  have  large  variances  and 
bias  values  when  completed.  Having  a  priori  knowledge  of  the  maximum  possible  area  (0) 
the  feature  can  move  in  during  a  sample  interval,  eliminates  the  need  for  nonparametric 


methods.  The  problem  simplifies  to  one  of  linear  mapping  and  classification.  Dropping  the 
a  posteriori  term,  subsequent  searching  is  dependent  on  a  Maximum  Likelihood  Estimation 
(MLE): 


where  A,  is  the  search  area  in  pixels,  and  Rvmax  is  the  maximum  distance  a  feature  could 
travel  in  one  direction  given  its  maximum  expected  velocity  and  the  sampling  period. 
Rvmax  is  determined  by; 

,  (2.8) 

where  Vmax  is  the  maximum  expected  velocity  of  the  feature  in  meters  per  second,  t  is  the 
time  in  seconds,  and  Ip  is  the  side  dimension  of  the  pixel  in  meters. 


Figure  2.6:  Ideal  feature  space  mapping  from  pattern  tile  •  search  tile  comparison. 
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When  the  search  area  is  sampled,  patterns  will  map  into  two  feature  spaces:  a 
space  where  the  minimum  criteria  are  met,  wi,  and  a  space  they  do  not,  u;2  (Figure  2.6). 
This  sorting  is  accomplished  using  classifiers  (Figure  2.7)  to  discriminate  between  groups. 
The  criteria  used  in  this  work  are  the  statistical  correlation  of  the  pattern  tile  -  search  tile 
pair,  limiting  the  root  mean  square  (RMS)  temperature  variation  between  the  search  and 
pattern  tiles,  and  limiting  the  pattern  and  search  tile  inter-temperature  ranges.  Selection 
of  the  matching  feature  is  made  from  the  ufi  space,  picking  the  feature  which  best  matches 
the  pattern  tile. 


Figure  2.7:  Search  area  pattern  sorting  and  selection  using  classifiers. 

Selection  of  the  best  match  is  a  forward  selection  process  of  the  members  of  the  u>i 
space  using  a  basic  branch  and  bound  procedure.  The  number  of  possible  matches  in  a 
searched  area,  M,  is  defined  by: 


M  =  As- 1. 


(2.9) 


The  number  of  features  which  meet  the  classification  criteria,  M,  is  determined  by  the 
number  of  features  discarded  during  classification,  ft,  by: 


M  =  M  -  n  . 


(2.10) 
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Final  selection  of  the  optimum  feature  Mopt  assumes  monotonicity  of  the  samples,  which 
allows  the  process  to  majumize  the  samples  and  pick  the  best  feature  such  that: 

Mopt  >  Mj»-i  >  >  Ml .  (2.11) 


Maximum  ocean  surface  current  flow  rates  are  known  for  most  of  the  world,  seldom 
exceeding  one  meter  per  second  in  the  open  ocean,  and  less  than  two  meters  per  second 
for  the  Delaware  Bay  and  estuary.  Using  the  maximum  likelihood  -  linear  classiflcation 
approach,  search  areas  are  limited  to  the  region  around  each  point  which  is  most  likely  to 
contain  the  matching  vector. 

2.2.1. 1  Search  area  dimension  determination 

The  search  area  dimension  is  treated  as  a  variable  in  this  study.  It  is  calculated  for 
each  image  pair  processed,  depending  on  the  following: 

•  Image  pixel  dimension. 

•  Time  elapsed  between  images. 

•  Expected  flow  rate  in  the  search  area. 

The  image  pair  date  time  group  is  used  to  calculate  elapsed  time  between  images  in 
seconds.  This  time  is  multiplied  by  the  expected  flow  rate  (meters  per  second)  and  the 
product  is  divided  by  the  pixel  dimension  (1100  meters  for  this  study)  to  generate  the 
search  radius  vector.  A  notional  square  is  calculated,  using  twice  the  search  radius  vector 
as  its  side  dimension.  This  box  is  centered  on  the  geographic  center  of  the  pattern  tile, 
and  the  search  is  started  (Figure  2.8). 

2.2.2  Feature  recognition  and  selection 

2. 2.2.1  Feature  recognition 

Search  tiles  are  matched  to  pattern  tiles  using  a  common  statistical  correlation 
algorithm  which  has  been  constrained  in  a  manner  similar  to  SEDA  algorithm.  StatisticaJ 


Figure  2.8:  Search  areas  and  best  match  selection,  continental  shelf  region,  0132  hours 
local,  30  June  1989.  Pattern  tile  positions  are  marked  in  blue,  optimum 
match  tiles  are  black. 

correlation  treats  the  pattern  and  search  tiles  as  real,  discrete,  two  -  dimensional  func¬ 
tions,  p(m,  Ti)  and  s(m,  n).  From  Wahl  and  Simpson  [3],  the  autocovariance  and  cross  - 
covariance  are  given  by: 

C„  =  n)  -  77,}{5(m  -f-  mo,  n  +  no)  -  t/,}]  ,  (2.12) 

C,p( mo,  no)  =  E[{s(m,n)~T),}{p{m+mo,n  +  no)-r]p}],  (2.13) 
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where  £?[•]  is  the  expected  vadue,  (mo,  no)  is  the  spatial  lag  between  functions,  rf,  and  rjp 
are  given  by  the  function  means: 


Tit  -  i;[s(m,  n)] , 

(2.14) 

rip  =  i;[p(m,n)]. 

(2.15) 

The  function  variances  and  Cp  in  fact  are  the  zero  -  lag  autocovariances.  This 
allows  the  correlation  coefficient  to  be  defined  as: 


r,p(mo,  no) 


Ctpjmo,  no) 

(T,(Tp 


(2.16) 


such  that  ||r«p(mo,  no)||  <  1  and  rM(0,0)  =  1.  This  means  that  if  a  chosen  search  tile  is 
an  exact  spatially  dislocated  version  of  the  first  signal,  then: 

s(m,  n)  =  p(m  +  mo,  n  +  no) ,  (2.17) 


and  r,p{mo,  no)  must  equal  one.  In  matrix  form,  the  correlation  coefficient  can  be  calcu¬ 
lated  using  the  following: 


r,pik,l)  = 


Ei  EjM*  +  k,j  +  l)-  r),{k,  l)]\p{i,j)  -  Vp] 

[a  E, ■(»(•'  +  k,j  +  i)-  01=  EvEilKi, ;')  -  ^ 


(2.18) 


Applying  tile  correlation  in  this  manner  estimates  the  image  spectrum  in  an  au¬ 
toregressive  manner  [10].  The  advantage  in  using  this  method  as  a  selection  criteria  is  that 
it  is  less  sensitive  to  background  noise.  Use  of  the  correlation  method  to  select  matching 
features  is  well  accepted  and  is  widely  used.  While  computationally  intensive,  for  pattern 
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tiles  12  X  12  and  smaller,  it  has  been  proven  that  calculating  the  correlation  directly  is 
faster  than  using  either  discrete  or  fast  fourier  transformation  techniques  [3]. 

2. 2.2.2  Feature  rotation 

Given  time  periods  between  images  of  12  -  24  hours,  it  is  almost  certain  that  ocean 
thermal  features  are  rotating  as  well  as  displacing  in  the  X  -  Y  plane.  This  is  especially 
true  in  eddy  flow  areas  such  as  the  one  observed  over  the  continental  shelf,  south  of  Long 
Island.  Methods  exist  [11]  to  transform  the  pattern  tile  to  different  axis  orientations. 
A  single  rotation  transformation  to  allow  for  a  left  or  right  movement  would  increase 
computations  by  a  factor  of  three.  Rotation  transformations  are  not  calculated  in  this 
program  to  speed  execution  time. 

Emery  [12]  states  that  a  correlation  above  0.4  is  has  a  99%  significance  for  features 
that  have  rotated  in  the  X  -  Y  plane  as  much  as  40  degrees.  Wahl  and  Simpson  [3] 
recommend  correlations  of  0.65  or  greater.  This  study  uses  a  minimum  coefficient  of 

T,p  =  0.60. 


2.2. 2.3  Feature  selection 

To  reduce  the  number  of  calculations  performed,  a  tile  with  dimensions  equal  to 
the  edge  detection  tile  is  first  moved  through  the  search  area.  A  correlation  between  this 
tile  and  the  center  of  the  pattern  tile  is  made  and  compared  to  a  threshold.  If  the  desired 
value  is  met,  the  full  search  tile  area  is  remapped  around  the  center  of  the  smaller  tile. 

The  search  tile  area  is  then  correlated  with  the  full  pattern  tile.  If  the  correlation 
criteria  is  met,  the  root  mean  square  (RMS)  difference  of  the  tiles  is  evaluated  to  insure  that 
the  mean  temperature  change  of  the  tile  does  not  exceed  a  physically  realistic  maximum. 
This  limit  is  typically  in  the  range  of  3  -  5°  C  (24  -  40  pixel  counts).  Additionally,  the 
highest  and  lowest  ranked  pixels  are  also  compared  to  the  temperature  thresholds. 

These  error  checks  are  required  because  the  correlation  coefficient  is  affected  by 
spatial  dislocation  in  three  dimensions,  not  just  two  (Figures  2.9,  2.10  and  2.11  ).  A 
search  tile  that  is  an  exact  match  for  the  pattern  tile  in  the  X  and  Y  dimensions,  can 
also  vary  in  the  Z,  or  temperature  dimension,  as  well.  A  uniform  change  in  temperature 
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Figure  2.9:  Pattern  tile  -  search  tile  correlation:  r,p  =  0.85. 

can  yield  a  high  correlation,  but  give  an  erroneous  selection.  This  occurs  most  commonly 
over  land  and  water  in  the  presence  of  light  or  scattered  clouds.  Light  clouds  may  differ 
in  signal  by  30  pixel  counts  or  less  from  a  water  surface. 

As  the  program  finds  features  that  meet  the  threshold  criteria,  they  are  compared 
to  the  currently  stored  value.  If  the  feature  has  a  higher  correlation,  it  overwrites  the 
previous  entry  and  increments  a  counter.  When  the  search  area  is  exhausted,  the  area 
with  the  highest  correlation  to  the  pattern  tile  is  used  to  compute  the  surface  flow. 


26 


Figure  2.10:  Pattern  tile  -  search  tile  correlation:  r^p  =  0.63. 

2.3  Direction  and  magnitude  computation 

Once  search  tile  selections  have  been  made,  the  x  and  y  values  are  compared  for 
each  tile  pair.  Trigonometric  calculations  are  performed  to  determine  vector  length  (mag¬ 
nitude),  direction  (azimuth)  and  flow  rate  in  meters  per  second.  To  improve  computational 
speed,  the  calculations  are  switched  to  specific  routines  based  on  the  difference  of  the  tile 
locations. 

The  magnitude  of  the  vector  is  computed  the  expression: 

1  =  y/ixo-  ®i)^  +  (yo  -  yiy .  (2.19) 

The  direction  is  found  using: 
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Figure  2.11:  Pattern  tile  -  search  tile  correlation:  r,p  =  0.25. 


-  hsum 

6  =  arccos  — ; — 


(2.20) 


ff  is  then  converted  to  azimuth  based  on  the  quadrant  in  which  the  angle  was  computed. 

Trigonometric  transformations  for  arrowheads  are  also  made  in  these  subroutines, 
using  the  expression: 


^2(arrow)  =  ±  COS X  M)  ±  ±  sin^^  X  M)  ,  (2.21) 

where  is  the  angle  of  the  arrow  head  from  the  line,  and  M  is  the  arrowhead  length. 
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The  vector  lines  and  arrowheads  are  drawn  on  the  image  using  a  C-language  version 
of  Bresssenham’s  integer  line  drawing  algorithm  [13]. 

2.4  Cloud  detection 

The  ability  of  the  SEDA  and  constrained  correlation  methods  to  operate  with 
images  that  contain  land  and  clouds  is  a  significant  advantage  (Figure  2.12).  Of  more 
than  131  images  reviewed  from  the  period  1  March  to  30  June  1989,  only  6  sets  were 
sufficiently  open  in  the  area  of  interest  to  permit  study. 


Figure  2.12:  Pattern  tile  selection  in  cloudy  conditions,  continental  shelf  region.  1308 
hours  local,  29June  1989. 

However,  the  SED.A.  method  is  not  adequate  as  a  method  for  screening  images  in 
an  automated  manner.  Uniformly  thin  cloud  layers  can  defeat  the  method,  producing  a 
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laxge  number  of  false  vectors.  To  screen  large  numbers  of  images,  a  more  effective  cloud 
screening  method  such  as  Simpson’s  [14]  or  Gautier’s  [15]  is  recommended. 


Chapter  3 


SURFACE  FLOW  IN  THE  DELAWARE  COASTAL  REGION 

3.1  Area  studied 

This  project  observed  the  surface  thermal  feature  displacement  in  two  areas:  the 
Delaware  coastal  region,  and  an  area  of  interesting  eddy  flow  over  the  continental  shelf, 
south  of  Long  Island.  Comparison  to  in  situ  anchored  buoy  data  has  been  made  for  the 
Delaware  coastal  region. 

Images  used  for  this  study  are  listed  in  table  3.1  by  their  date  time  group  (DTG). 
The  DTG  is  can  be  read  as  follows: 

Example  DTG:  89181063256 
yy  ddd  hh  mm  ss 

where  yy  is  the  year  (89  =  1989),  ddd  is  the  Julian  date,  hh  is  the  hour  in  Greenwich  Mean 
Time  (GMT),  mm  are  the  minutes  and  ss  are  the  seconds. 

3.2  Analysis  of  surface  flow  in  the  Delaware  coastal  region 

The  flow  of  surface,  subsurface  and  bottom  currents  in  the  Delaware  coastal  region 
has  been  extensively  studied.  A  recent  report  by  Miinchow  [16],  provides  an  excellent 
presentation  of  this  information.  Figure  3.1  shows  approximate  locations  for  21  of  31 
anchored  buoys  used  in  his  study.  These  buoys  have  mean  observations  for  current  velocity 
and  direction  for  depths  between  5  and  10  meters.  Satellite  observations  in  the  coastal 
area  have  been  correlated  with  these  observations  as  listed  in  Table  3.2. 
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Table  3.1:  NOAA  satellite  images  used  in  this  study,  listed  by  date  time  group. 


NOAA  Images  Used 

Period 

Image  DTG 

10-11  March  1989 

89069170021 

89069184155 

89070183134 

28-29  May  1989 

89148071223 

89148165803 

89149070111 

10-12  June  1989 

89161063805 

89161180325 

89162062806 

89163061744 

89163080102 

29-30  June  1989 

89180180847 

89181063256 

89181175808 

3.2.1  Observations  10  March  >  30  June  1080 

3.2.1.1  10-11  March 

Figure  3.2  shows  that  thermal  features  are  detected  in  the  coastal  region.  The 
program  performs  well  in  selecting  only  water  features.  However,  no  search  tiles  were 
selected  for  either  image  pair  from  this  period. 

The  output  data  shows  that  the  sea  surface  temperature  (SST)  and  bay  temperature 
is  very  low,  5.25®  to  5.62°  C  respectively.  This  represents  a  range  from  26  to  61  in  pixel 
values.  Of  the  31  pattern  tiles  selected,  16  found  no  reduced  set  search  tiles  that  met  the 
correlation  criteria.  The  remaining  15  computed  from  79  to  137  full  tiles,  but  none  had 
an  Ttp  >  .50. 

The  low  correlation  suggests  a  very  deep  mixed  layer  due  to  the  cold  water  tem¬ 
perature  at  the  surface.  This  means  that  the  vertical  diffusion  condition  (fj  <  I7)  is 
not  met  due  to  continuous  vertical  mixing  of  surface  and  deep  (greater  than  50  meters) 
water  [17]. 
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Table  3.2:  Selected  anchored  buoy  near  surface  current  measurements  (from  Miinchow). 


M2  tidal  properties 


Mooring 

Depth 

(meters) 

Axis 

degrees 

Velocity 

cm/s 

1 

4 

127 

94.3  ±  0.3 

2 

4 

123 

80.2  ±0.1 

3 

6 

126 

60.8  ±  0.2 

4 

5 

139 

43.5  ±0.1 

5 

5 

137 

69.6  ±0.1 

6 

5 

120 

24.1  ±  0.4 

7 

5 

148 

21.0  ±0.1 

8 

5 

160 

18.7  ±0.3 

9 

5 

160 

16.5  ±  0.2 

11 

6 

106 

17.6  ±  0.2 

12 

3 

129 

24.2  ±  0.2 

14 

6 

152 

25.1  ±  0.0 

16 

6 

97 

17.2  ±0.3 

17 

6 

97 

16.8  ±  0.4 

18 

6 

111 

19.9  ±0.3 

19 

6 

202 

7.1  ±0.3 

20 

6 

187 

8.2  ±  0.4 

21 

6 

179 

10.8  ±0.4 

22 

8 

150 

21.6  ±0.2 

24 

7 

157 

7.7  ±0.2 

25 

8 

144 

9.4  ±  0.3 
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Figure  3.1:  Selected  anchored  buoy  locations  for  the  Delaware  coastal  region,  which 
were  used  to  evaluate  accuracy  of  satellite  derived  surface  current  flow. 

3.2. 1.2  28-29  May  1989 

Figure  3.3  shows  magnitude  and  direction  vectors  for  a  12  hour,  3  minute  period, 
28  May  1989.  The  three  images  available  were  processed  as  two  pairs.  Only  the  pair 
89148165803  to  89149070111  were  acceptably  cloud  free  to  observe  the  Delaware  coastal 
region,  producing  22  pattern  tiles  and  10  vectors. 

The  continental  shelf  area  was  observed  using  the  grid  technique  for  both  days.  In 
each  case,  49  pattern  tiles  were  selected,  with  only  10  search  tile  selections  for  the  period 
from  07-1600  GMT,  and  16  selections  for  the  16-0700  GMT  period.  The  small  selection 
rate  is  attributed  to  a  combination  of  vertical  mixing  due  to  low  water  temperature  (15.5° 
C  in  the  bay,  12.6°  C  offshore)  and  rotation  of  features  in  the  eddy  areas. 


Figure  3.2:  Pattern  tile  selection  in  the  Delaware  Bay  and  coastal  region.  1100  hours 
local,  10  March  1989, 

3.2. 1.3  10  -12  June  1989 

Figure  3.4  shows  a  high  density  sampling  of  the  flow  around  the  Delaware  coastal 
region,  12  June  1989.  Four  images  were  available  for  analysis,  including  one  pair  with  a  24 
hour  spacing  and  another  pair  with  a  2  hour  interval.  In  four  comparisons,  94  vectors  were 
produced  from  a  pattern  tile  set  of  142.  This  increase  in  selection  shows  stratification  has 
occurred  and  a  surface  mi.xed  layer  has  formed  and  is  more  stable  with  water  temperatures 
reaching  18.63°  C  in  the  bay  and  16.13°  C  offshore. 

The  continental  shelf  area  was  evaluated  using  both  grid  and  SED.\  methods.  93 
out  of  187  pattern  tiles  selected  produced  vectors.  The  flow  pattern  was  not  well  defined 
in  this  set,  but  is  similar  to  the  flow  defined  in  Figure  3.6. 


Figure  3.3:  Labeled  flow  vectors  in  the  Delaware  Bay  and  coastal  region  for  a  12  hour 
period  ending  0201  hours,  28  May  1989. 
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Figure  3.4:  High  density  flow  vectors  in  the  Delaware  coastal  region,  for  a  2  hour,  16 
minute  period  ending  0201  hours  local,  12  June  1989. 


Figure  3.5;  Surface  flow  vectors  in  the  Delaware  coastal  region,  for  an  11  hour,  2S 
minute  period  ending  1258  hours  local,  30  June  1989. 

3.2.1.4  29-30  June  1989 

Figure  3.5  shows  general  surface  flow  in  the  Delaware  coastal  region,  30  June  1989. 
Three  images  were  available  for  this  period.  Each  was  spaced  by  appro.\imately  12  hours, 
and  were  processed  cis  two  pairs.  105  vectors  were  obtained  from  171  pattern  tiles.  Search 
tile  matching  was  hampered  by  light  clouds  which  moved  into  the  offshore  region  in  each 
comparison  image.  Sea  Surface  temperature  was  very  warm,  measured  at  21.38°  C  for  the 
bay,  and  20.5°  C  for  the  offshore  area. 

The  continental  shelf  region  was  evaluated  using  both  the  SEDA  and  grid  methods. 
Figure  3.6  shows  .53  vectors  derived  from  110  pattern  tiles.  Feature  rotation  is  the  most 
probable  cause  for  the  low  selection  rate.  Horizontal  and  vertical  diffusion  are  certainly  a 


Figure  3.6:  Surface  flow  vectors  in  the  continental  shelf  region,  for  an  11  hour,  26  minute 
period  ending  1258  hours  local,  30  June  1989. 

factor  in  the  eddy  regions. 

3.2.2  Comparison  of  satellite  observations  and  buoy  data 

Verification  of  this  program  to  accurately  measure  displacement  of  surface  thermal 
features,  and  judge  the  validity  of  using  this  movement  to  infer  actual  near  surface  currents, 
was  accomplished  by  comparing  program  derived  vectors  with  in  situ  buoy  data.  110 
vectors  were  matched  to  17  buoy  positions.  The  mean  of  these  values  for  each  station  was 
derived  (Table  3.3),  and  the  results  were  correlated  with  good  results  (Table  3.4). 

These  values  were  obtained  by  computing  the  mean  of  vectors  obtained  by  visually 
matching  program  vectors  with  the  closest  buoy  station.  Vectors  were  rejected  if  the 


39 


Table  3.3:  Mean  satellite  observations  for  selected  m  situ  buoys  from  Table  3.2  Val¬ 
ues  were  determined  by  matching  satellite  derived  vector  origins  with  buoy 
positions.  The  values  listed  are  the  means  of  multiple  observations  for  each 
location. 


Mean  satellite  -  buoy  observations 

Mooring 

Axis 

degrees 

Velocity 

cm/s 

1 

131.02 

112.80 

4 

134.37 

36.75 

5 

159.86 

79.00 

6 

116.45 

34.00 

7 

90.00 

27.0 

9 

174.34 

12.50 

11 

108.43 

42.00 

12 

165.85 

26.00 

14 

180.19 

32.33 

16 

92.33 

26.00 

17 

111.07 

27.00 

18 

117.96 

23.67 

19 

217.46 

12.50 

20 

193.28 

18.50 

21 

168.74 

16.40 

22 

146.48 

40.00 

24 

164.52 

10.50 

Table  3.4:  Satellite  -  Buoy  azimuth  and  velocity  comparison. 


Mean  satellite  - 

buoy  comparison 

Parameter 

Result 

Azimuth  Correlation 

.83 

Std.  Deviation  {<tax) 

20.45“’ 

Mean  Error 

-4.38° 

Velocity  Correlation 

.96 

Std.  Deviation  {av) 

7.76  cm/s 

Mean  Error 

-7.75  cm/s 
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magnitude  or  direction  were  obviously  false  (velocity  errors  greater  than  a  factor  of  3,  or 
more  than  60°  oif  in  azimuth),  due  to  cloud  or  land  contamination  (see  vector  11,  Figure 
3.3),  prior  to  determining  the  mean  value.  Vectors  with  direction  errors  of  180°  ±  60°  were 
considered  to  be  on  the  flow  axis  and  were  converted  by  adding  or  subtracting  180°  as 
appropriate.  Of  the  110  vectors  considered,  27  were  rejected  using  this  criteria,  giving  a 
final  sample  size  of  83. 

The  agreement  in  velocity  between  the  satellite  and  buoy  measurements  is  partic¬ 
ularly  good  with  a  av  of  only  7.76  cm/s  over  a  range  from  9  to  112  cm/s.  The  azimuth 
correlation  is  also  a  good  match.  The  of  20.45°  reflects  the  smaU  sample  size,  and 
would  decrease  with  more  data  points. 

3.2.3  Comparison  with  the  Simpson  Minimum  Distortion  Log  Search  method 

A  program  using  the  Simpson  Minimum  Distortion  Log  Search  Method  [6]  is  avail¬ 
able  on  line  in  the  University  of  Delaware  Center  for  Remote  Sensing.  In  its  current  form 
only  the  two  images  from  the  30  June  1981  period  are  available  for  comparison,  using  a 
40  X  40  pixel  grid  in  the  continental  shelf  region.  This  program  was  run  in  the  grid  mode 
producing  13  vectors  from  49  pattern  tiles  selected.  The  correlations  for  these  methods  is 
very  low:  .31  in  velocity,  and  .23  in  azimuth. 

The  most  probable  reasons  for  these  low  correlations  are  a  differing  of  search  area 
dimensions  (this  program  searches  an  area  four  times  larger  than  the  Simpson  method), 
and  the  difference  in  search  tile  selection  methods.  The  Simpson  program  forces  selection 
of  a  vector  for  each  pattern  tile  in  the  area  which  shows  the  least  distortion.  This  program 
is  more  stringent,  requiring  a  statistical  correlation,  rejecting  features  which  do  not  meet 
desired  criteria. 


Chapter  4 


CONCLUSION 


4.1  Summary 

This  program  is  a  tool  which  can  be  used  to  estimate  surface  and  near  surface  cur¬ 
rents,  by  objectively  measuring  the  displacement  of  large  water  surface  thermal  features. 
It  utilizes  commonly  available  satellite  images  with  no  special  enhancement  or  processing 
requirements.  The  constrained  Ordered  Statistical  Edge  Detection  algorithm  employed 
to  select  features  is  capable  of  discriminating  between  desired  water  surface  features  and 
other  background  objects  such  clouds  or  land.  The  program  is  written  in  the  C  computer 
programming  language,  and  has  been  run  on  a  variety  of  common  systems  including: 

•  Silicon  Graphics  -  Personal  Iris. 

•  Digital  Equipment  Corporation  -  DEC  Station  5000. 

•  Sun  Microsystems  -  Sun-3  operating  system. 

The  program  produces  output  in  image  and  text  form.  Images  are  labeled  with  identifi¬ 
cation  numbers  to  allow  comparison  of  the  text  file  and  the  output  images. 

The  high  correlation  values  derived  from  comparing  satellite  image  derived  mea¬ 
surements  with  in  situ  buoy  data,  combined  with  the  ability  to  use  images  which  have 
clouds  and  land  present,  make  this  program  a  good  tool  for  the  study  of  coastal  and  open 
ocean  waters.  Results  generated  by  this  program  will  provide  useful  input  for  mixed  layer 
depth  modeling  [18]. 
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Appendix  A 


COMPUTER  PROGRAM  DESIGN 

The  computer  program  which  generated  the  data  and  images  used  in  this  work  has 
an  operating  system  consisting  of  seven  major  sub-groupings: 

•  Program  Setup  and  User  Interface. 

•  Image  Input. 

•  Pattern  Tile  selection. 

•  Search  Tile  /  Pattern  Tile  matching. 

•  Surface  Flow  Computation. 

•  Image  Graphics. 

•  Output. 

Program  operation  is  illustrated  in  figure  A.l. 

A.l  Progrsun  setup  and  user  interface 

Program  ^itialization,  setup  and  user  interface  are  provided  by  the  following  rou¬ 
tines: 

•  main.c. 

•  setup.c. 

•  userinput.c. 
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^ - Multiple  Image  Processing 

'*■  ►  Program  /  User  interaction 

.  Program  status  updates 


Figure  A.l:  Ocean  flow  estimation  program  flow  chart. 

•  date.c. 

•  characters  .c. 

•  image.h. 


Main.c  is  a  standard  C  routine  which  defines  the  computer  program  fiow  chart. 
Setup.c  contains  the  default  values  for  the  program  global  variables.  Userinput.c  pro¬ 
duces  a  series  of  screen  menus  and  prompts  to  allow  single  run  changes  to  the  program 
parameters.  Date.c  strips  the  date  time  groups  from  the  command  line  starting  the  pro¬ 
gram  to  compute  elapsed  time.  Characters.c  contains  the  character  bit  map  for  graphic 
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labels  placed  on  the  processed  images.  Image.h  is  a  header  file  which  contains  structure 
variable  and  global  function  and  value  declarations. 

A. 2  Input  /  Output  functions 

Image  input  and  output  for  this  program  are  fairly  standard  routines  that  are  in 
common  use.  The  input  functions  are  found  in; 

•  get  image.c. 

•  get  image2.c. 

Input  images  may  be  processed  without  headers  or  in  the  Interactive  Viewing  and 
Analysis  System  (IVAS)  [19]  format,  concatenated  with  a  four  line  header. 

Output  images  can  be  output  with  or  with  the  header  by  changing  commenting 
around  the  header  write  commands  in  the  output  files.  Headers  must  be  removed  when 
using  Ximage  or  similar  non-IVAS  compatible  viewing  tools.  Image  output  is  performed 
in  these  files: 

•  put  image.c. 

•  put  image2.c. 

•  put  image3.c. 

•  put  image4.c. 

•  put  imageS.c. 

Text  files  are  generated  at  the  end  of  computations  from  the  routines: 

•  data.c. 

•  matlab  out.c. 

This  output  can  be  customized  to  meet  user  needs.  Current  output  goes  to  two 

files: 
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•  flowdata[n]. 

•  matdata[ii]. 


flowdataln]  is  a  text  file  that  lists  results  of  each  run  (see  Appendix  D).  The  counter 
[n]  is  incremented  when  multiple  images  are  processed  to  provide  the  data  from  each  image 
pair  comparison.  matdata[n]  is  a  matrix  output  of  the  pattern  tile  -  search  tile  pairs  in 
PRO  -  MATLAB  format  [20]. 

Userinput.c  is  a  routine  which  allows  the  user  to  make  one  time  changes  to  the 
program  parameters.  This  routine  also  prompts  the  user  for  setup  information  which  may 
frequently  change.  Shown  below  is  an  example  of  the  user  screen  if  all  options  are  reviewed 
and  changed. 


>  flow  1  80089170021  8906918416B 

Are  the  laaga  diaenaioni  (Xdia)  S12  X  612  (Ydim)?  (y/u)  :y 
Do  you  vant  to  vlov  dolault  setup  ?  (y/u):  y 

Changes  are  sMide  lor  this  session  only,  changes  in 

default  values  are  aade  by  editing  ’setup. c' 

Area  of  Interest:  ZBin-40  ynin-210  ,  xmx-160  ynaz-300 

Pattern  tile  Default  selection  criteria: 

Search  Option: 

Tile  size  : 

Edge  detect  tile  size: 

Ptile  Biniaran  offset: 

Filter  aininuB  edge  threshold: 

Pixel  difference  value  range: 

Search  tile  Default  selection  criteria: 

Search  area  flow  rate:  0.60 

Stile  search  increaent:  1 

Allowed  velocity  error:  0.20 

Allowed  Tenperature  error:  4.00 
Allowed  variance  froa  aean:2.00 
Niniana  correlation  value:  0.60 

Changes  can  be  aade  by  entering: 

a  -(area  of  interest)  p  -(ptile)  s  -(stile)  d  -(done):a 


edge 

11 

6 

20 

6 

12 
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Eater  search  area  rain:  40 
Enter  search  area  ynin:  210 
Enter  search  area  max:  160 
Enter  search  area  ymaz:  300 

Enter  a  to  re-edit,  p  -ptile,  s  -stile,  or  d  -done:p 

Enter  Search  option  (l^Edge  Detection,  2sGrid  Search):  1 

Eater  Ptile/Stile  size  number(Larger  than  Edge  k  odd;  6,7,9,11):  11 

Enter  Edge  Detect  tile  size  (should  toe  odd;  3,6,7...):  6 

Enter  Ptile  ainimun  offset:  20 

Enter  filter  sdninna  edge  threshold:  6 

Enter  pixel  difference  value  range:  12 

Enter  p  to  re-edit,  a  -area  of  int,  s  -stile,  or  d  -done:8 

Enter  search  area  flow  rata:  .6 
Enter  Stile  search  increment:  1 

Enter  allowed  Tile  RMS  temp,  difference  (  deg.  C):  4 
Enter  minimum  correlation  value  (0.1  -  1.0,  usually  0.6):  .6 

Enter  s  to  re-adit,  a  -area  of  int  p  -ptile,  or  d  -done:d 


Search  Area  «  4,  for  a  max  flow  rate  of  0.6 

Do  the  input  images  have  headers?  (y/n):n 

Do  you  want  the  output  images  to  have  headers?  (y/n):a 

Min  mask  value:  26,  Max  mask  value:  61,  Del  bay:  46  (6  C),  Coast:  42  (6  C) 

Ptiledata  start 

ptiledata  complete  10  ptiles 

starting  stile. c 

finished  with  stile. c 

writing  boxes  onto  st  image 

writing  stiles  to  image 

Starting  vector. c 

Output  is  going  to  flowdatal 

Output  is  going  to  natdatal.a 

> 


A.3  Pattern  tile  selection 

Pattern  tiles  are  selected  in  the  iiie: 


e  fUter.c  (ordered  statistic  edge  detection). 
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•  grid.c  (Specified  x  -  y  placement). 

Pattern  tile  dimensions  are  variable.  For  this  study  a  5  x  5  array  (  25  subtiles  )  was 
used  to  detect  features,  an  11  x  11  tile  mapped  features  that  were  selected.  The  values 
for  the  feature  discrimination  values  maski  and  mask2  are  calculated  in  the  routine: 

•  maskvalue.c 

Masks  can  be  placed  anywhere  on  an  image  depending  on  user  needs.  For  this 
study,  masks  were  placed  in  the  Delaware  bay  and  40  kilometers  off  shore  from  Cape  May, 
New  Jersey.  The  masks  are  computed  from  tiles  the  same  size  as  the  mapped  pattern  tile. 

A.4  Search  tile  operations 

Search  tile  correlation  and  selection  take  place  in  the  file: 

•  stile.c. 

A.4.1  Correlation  coefficient  and  RMS  difference  computation 

To  reduce  the  amount  of  dynamic  memory  required,  only  search  tiles  that  meet  the 
classification  criteria  are  stored  temporarily  in  a  two  dimensional  structure.  If  a  search 
tile  fails  to  meet  correlation  or  RMS  difference  criteria,  it  is  not  stored  and  the  count 
is  not  advanced.  Additionally  if  the  corrdation  value  is  not  greater  than  the  previously 
recorded  matched  tile,  the  tile  is  not  recorded  and  the  count  is  not  advanced.  This  means 
that  when  the  search  tile  box  has  been  exhausted,  the  closest  match  is  at  the  top  of  the 
structure  and  is  recorded  for  vector  calculation. 

Using  a  reduced  set  search  (5x5  versus  11  x  11  )  reduces  the  total  number  of 
calculations  by  roughly  a  factor  of  four.  A  30  x  30  search  box  using  a  5  x  5  search  tile 
has  625  iterations  of  25  calculations  each  (15,625  calculations).  Searching  with  an  11  x  11 
tile  would  take  75,625  csilculations.  The  actual  total  lies  somewhere  in  between.  Table 
A.l,  shows  numbers  from  a  typical  run.  The  small  number  of  times  that  the  full  tile  was 
computed  produced  a  factor  of  four  improvement  in  execution  speed. 

If  no  search  tile  meets  either  of  the  threshold  criteria,  the  pattern  tile  location 
values  are  recorded  for  the  search  tile  and  no  fiow  magnitude  or  direction  is  calculated. 
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Table  A.l:  Reduced  set  search  calculation  savings 

Search  Box  Dimension  47  x  47 
_ 2209  iterations  per  tile _ 

Tiles  Full  Tiles  Cost  Cost  Savings 
Selected  Computed  all  11  x  11  this  run  this  run 
15  473  4,009,335  940,833  76.54  % 


A.5  Vector  magnitude  and  direction  calculation 
Vector  operations  are  performed  in  the  file: 

•  vector.c. 

Once  search  tile  selections  have  been  made,  the  x  and  y  values  are  compared  for  each 
p-tile  -  8-tile  pair.  Trigonometric  calculations  are  performed  to  determine  vector  length 
(magnitude),  direction  (azimuth)  and  flow  rate  in  meters  per  second.  To  improve  compu¬ 
tational  speed,  the  calculations  are  switched  to  specific  routines  based  on  the  difference  of 
the  p-tile  and  s-tile  locations: 

hsum  =  s  -  tile{x)  —  p  -  tile{x) , 


vsum  =  s  —  tile{y)  —  p  —  tile{y) . 
Magnitude  of  the  vector  is  computed  the  expression: 

I  =  Vhsum?  -I-  /sum*  . 


Direction  is  found  using: 


9  =  arccos 


hsum 

~1~ 


9  is  converted  to  azimuth  based  on  the  quadrant  for  with  the  angle  was  computed  (figure 
A.2). 
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Case4 


Case  1 


Case  2 


Figure  A.2:  Vector  transformation  diagram  used  to  calculate  vector  azimuth  and  end 
points  for  arrowheads 

Trigonometric  transformations  for  arrowheads  are  also  made  in  these  subroutines, 
using  the  expression: 

p2(arrow)  ~  (■^P2(iin*)  ^  cos^^  X  Af)  i  i  sin^  X  Af)  , 

where  0^  is  the  angle  of  the  arrow  head  from  the  line,  and  M  is  the  arrowhead 

length. 

The  vector  lines  and  arrowheads  are  drawn  on  the  image  using  a  C-language  version 
of  Bresssenham’s  integer  line  drawing  algorithm  [13]  which  is  found  in  the  routine: 

•  line.c. 

A.6  Program  validation 

To  check  the  program  to  insure  that  it  is  properly  mapping  and  calculating  pattern 
and  search  tiles,  the  program  was  run  using  the  same  image  for  tile  selection  and  compar¬ 
ison.  This  was  accomplished  by  copying  the  image  and  altering  its  date  time  group  by  12 
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hours.  The  program  was  executed  normally  with  the  desired  result  being  the  search  tile 
selection  of  the  pattern  tile.  Correlations  (r,p)  of  1.00,  changes  in  temperature  of  0.0,  and 
0.0  values  for  all  vector  calculations  indicate  a  working  program.  A  program  test  in  the 
Delaware  coastal  region  is  shown  below: 

From:  891810632S6.  To:  89181180000: 

Elapsed  time:  0  OayCs),  Total:  12  Hours,  32  Minutes,  56  Seconds 
Search  box  dim:61  Offset:  20  Tile  dim:  11  X  11  lumber  of  samples  9 
pt  Ingth(km)  az(deg)  flow  m/s 


0 

0.00 

0.00 

0.00 

1 

0.00 

0.00 

0.00 

2 

0.00 

0.00 

0.00 

3 

0.00 

0.00 

0.00 

4 

0.00 

0.00 

0.00 

5 

0.00 

0.00 

0.00 

6 

0.00 

0.00 

0.00 

7 

0.00 

0.00 

0.00 

8 

0.00 

0.00 

0.00 

Average  flou  distance:  0.00  km,  Avg  Az:  0.00  deg,  Avg  flow  rate:  0.00  m/s 


Search  tile  data  (z,y  coord  for  center  subtile) 


ptile  stile 

1  X 

1 

CC  val 

Pt  temp 

St  temp 

del  T 

Pt  sum  St  sum 

Ratio 

0 

116 

69 

241 

1.00 

22.65 

22.65 

0.00 

21826 

21826 

1.00 

1 

67 

129 

241 

1.00 

21.07 

21.07 

0.00 

20393 

20393 

1.00 

2 

146 

76 

262 

1.00 

20.71 

20.71 

0.00 

20044 

20044 

1.00 

3 

32 

109 

262 

1.00 

20.64 

20.64 

0.00 

19982 

19982 

1.00 

4 

29 

133 

262 

1.00 

21.67 

21.67 

0.00 

20877 

20877 

1.00 

6 

203 

91 

283 

1.00 

19.63 

19.63 

0.00 

18999 

18999 

1.00 

6 

26 

117 

283 

1.00 

21.68 

21.68 

0.00 

20989 

20989 

1.00 

7 

258 

90 

304 

1.00 

20.46 

20.46 

0.00 

19793 

19793 

1.00 

8 

30 

123 

304 

1.00 

22.22 

22.22 

0.00 

21613 

21613 

1.00 

System  working  values: 

Image  Xdim:  612  Ydim:  612,  Pizsize:  1100 

Search  Area  zmin:  40  ymin:  200  zmax:  140  ymax;  320,  Search  Option;  edge 

Ptile:  Edge  tile:  6,  Ptile:  11,  Offset:  20,  Grid  dim:  40,  Thresh:  6, 
Range:  12,  Naskl:  138,  Nask2:  196 
Areal  SST:  21.38,  Area2  SST:  20.25 

Stile:  search  dim:  40,  offset:  1,  Alloved  RMS  temp  change:  6.00 
Allowed  temp  range  from  tile  mean:  3.000,  CC  match(desired) :  0.600, 

Flow  rate(desired) :  0.600,  Max  allowed  flow:  0.800 
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Note  the  stile  column  in  Search  tile  data.  These  are  the  number  of  reduced  set  search 
tiles  (5  X  5)  which  met  the  selection  criteria  and  were  remapped  and  calculated  for  each 
pattern  tile. 

A.7  Image  graphics 

The  graphics  generated  by  this  program  covert  pixels  of  the  image  to  designated 
values.  This  makes  the  changes  permanent  as  far  as  the  output  images  are  concerned. 
This  is  done  to  facilitate  portability  of  the  output  images. 

Additionally,  a  small  pixel  character  map  was  written  for  this  program  to  draw  the 
date  time  group,  tile  labels  and  scales  shown  on  the  images. 

Currently  four  images  are  output  from  this  program: 

•  Pattern  Tile  Image  (P-tiles  written  on  the  input  image.  Each  tile  is  labeled  with  a 
text  file  identification  number). 

•  Search  Tile  Image  (  Comparison  image  overwritten  with  p-tiles,  selected  s-tiles,  and 
search  box). 

•  Vector  Image  (Comparison  image  overwritten  with  magnitude  and  direction  vectors 
,  Each  vector  labeled  with  text  file  identification  number). 

•  Arrow  Image  (Vector  Image  without  the  ID  numbers). 

Each  output  image  is  labeled  with  its  date  time  group  in  the  upper  left  corner, 
horizontal  and  vertical  0  -  30  kilometer  scale  in  the  lower  right  corner,  and  a  vector 
magnitude  arrow  in  meters  per  second. 

Files  which  perform  graphic  operations  are: 

•  graphics. c. 

•  label.c. 

•  line.c. 


•  veclabel.c. 
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A.8  Multiple  image  processing 

When  more  than  two  images  are  available,  they  may  be  sequentially  processed  to 
track  flow  over  time.  After  the  initial  two  images  are  evaluated  as  previously  described,  a 
routine: 

•  swapflle.c. 

is  called  which  moves  search  tile  data  into  pattern  tile  storage  locations  and  calls  the 
next  specifled  image  in  as  a  search  tile  image.  Search  areas  are  recomputed  based  on  the 
elapsed  time  between  the  new  image  pair,  and  the  pattern  search  and  vector  computations 
are  repeated.  At  the  completion  of  each  cycle,  text  and  record  images  are  generated. 


Appendix  B 


PROGRAM  CODE  EXTRACTS 


This  appendix  contains  extracts  from  the  current  operating  system  which  perform 
the  tasks  of  selection,  discrimination,  and  computation  of  surface  flow. 

B.l  Ordered  statistical  edge  detection  algorithm 

The  SEAD  algorithm  is  executed  in  two  parts:  a  reduced  set  search  to  locate  an 
acceptable  edge,  and  expanded  mapping  and  veriflcation  for  the  desired  pattern  matching 
dimension. 


B.1.1  Feature  detection  and  discrimination 

ka j+zdin*! ; 
conntsO; 

p8UB»0 ; 

lor(|sO;I<(edgedla) 

{ 

1  or  (M^O :  N<  (odgodia) ; 

i 

lsk‘«-l*zdiB-i-N; 

shellO [count] =iBage [1] ; 

ptilos [ptilocount] [count] . il  »  i; 

ptilos [ptilocount] [count] . j 1  =  j ; 

ptilos [ptilocount] [count] .kl  =  k; 

ptilos [ptilocount] [count] .sbolll  =  shellO [count] ; 

psuB=psuffl-*-shellO  [count] ; 

count+=l; 

> 

> 

ptilos [ptilocount] [0] .psual  =  psun; 

/*  storos  total  coll  value  for  match  in  stile. c  */ 

ptilos [ptilocount] [0] .avgvall  =  (int)((float)psuffl/(W)); 

/♦♦  rank  pts  ♦*/ 
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lor(q=0:q<(W-2) ;++q) 

f or(r=q+i ; r<W-l ; ++r) 

il ( shellO [q] >8hellO [r] ) 

■C 

temp=8hellOCq] ; 

8hellOCq]s8hellOCr] ; 

8li«110  Cr]  stamp; 

} 

} 

RangesshallO [W-l] -8hall0 [0] ; 

F6ature8  are  selected  here  lor  comparison. 

Thresh  is  the  threshold  r2aige  lor  selection 
as  an  edge.  A  value  ol  5  says,  that  when 
the  subtiles  are  ordered  by  value,  there  must 
be  at  least  a  value  ol  S  Irom  highest  to  lovest 

mg  is  the  maximum  separation  alloued  betveen 
min  and  max  subtile  values.  A  lou  value  but 
must  be  greater  than  Thresh,  discourages 
selection  ol  edges  ol  great  contrast  such  as 
eater  and  clouds. 

ma8kl,mask2  are  values  lor  areas  where  edges 
are  not  to  be  selected.  In  this  routine, 
the  value  ol  the  lowest  ordered  point 
must  be  higher  thm  the  min  mask  value  and 
the  max  value  must  be  lover  than  the  max  mask 

*•/ 

il(Range>Thre8h  At  Range<mg  A*  8hell0C0]>  maskl 
AA  8hellOCV]<ma8k2) 


/*  edge  detected  */ 
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B.1.2  Pattern  tile  mapping  and  verification 

Ksk-tilediff-tiledift*xdim: 

coimtsO : 

pstuBsO; 

for(l=0;l<(tiledim) ;++*) 

lor(M=0;M<(tiledim) ;++M) 

l=K-t-l*xdiB-t-N; 

shflllO [count] -image Cl] ; 

ptiles [ptilecount] [count] . i  =  i ; 

ptilesCptilecount] [count] .j  =  j; 

ptiles [ptilecount] [count] . k  =  k; 

ptiles [ptilecount] [count] . shell  =  shellO [count] ; 

psuB=psum-«-shellO  [count] ; 

count+Bl; 

} 

> 

ptilesCptilecount] [0] .psuB  =  psum; 

/*  stores  total  cell  value  for  match  in  stile. c  */ 
ptiles [ptilecount] [0] . avgval  =  (int)((lloat)psum/(WW)) ; 


/ee  Verify  that  all  points  in  mapped 
feature  lie  vithin  temperature 
range  **/ 

for(q*0;q<(WW) ;++q) 

f or (r=q+ 1 ; r<WW ; ++r) 

if ( ShellO [q] >shellO [r] ) 

temp=shellO[q] ; 

ShellO [q] ^shellO [r] ; 
shellO [r]=temp; 

> 

> 

if (shellO [0]>  maskl  kk  shellO [HW]<mask2) 
ptilecoiut-*--!; 

8S=1  ; 

/**  Apply  offset  to  sepaorate  points  **/ 
if  (j<jvalue) 

j=j+offset; 

} 
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else  if  (j>=j value) 
j=jBax; 

> 

> 


B.2  Subsequent  feature  recognition 

Feature  recognition  is  performed  by  performing  a  statistical  correlation  between 
reduced  set  tiles,  and  then  full  set  comparisons.  The  reduced  sets  are  compared  just  for 
correlation  and  mean  temperature  change.  K  they  meet  this  criteria,  the  full  pattern  tile 
dimensions  are  then  mapped  and  compared  for  correlation  and  points  are  ranked  to  insure 
that  no  points  exceed  the  allowable  temperature  range. 


B.2.1  Reduced  set  correlation 


k=j-t-i*zdin: 
countsQ ; 

88ua-0; 

for(|sO:I<(8dg«diB)  /*  88arching  each  8tile  */ 

i 

for(N>0;M<(edgedin) ;++M) 

l«:k-t-l*xdin-t-K; 

8helll [count] sinaga2 Cl] ; 

stile8 [a+l] [count] . 8helll  =  ahelll [count] ; 

stileaCa'fl]  [count]  .i  »  i* 

8tile8[a->-l]  [count]  .j  »  j; 
atileaCa+l] [count] .k  »  k; 

aaunsaaum+ahelll [count] ; 
count+*l ; 

> 

> 


/*  delining  valuea  lor  mean  avg  aum  match 
matchdels((w)e8um(ptilee8tile  -  8tUB[ptile]*8um[8tile])  / 
8qrt( (V) (aumCptile] *2-(8um[ptile] ) *2) )* 

( (W) (aumCatile] ~2-(8um[8tile] ) “2) )  */ 

aumysaaum; 

aumxsptilea [a] [0] . pauml ; 

nauaucysO; 

aumxaq^O; 

naumxaqsQ ; 
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•umsqxsO; 

ssaysqsO; 

nswiysq-O: 

suBsqysO; 

aunxysO; 


lor(c»0 ; c<W; ++c) 

{ 

Bimzy^snnxy+Catilas  Ca-t'13  Cc]  .sballlaptilasCa]  [c]  .ahelll) ; 

> 

n8UBXy=8unzy*(V} ; 

lor(c=0;c<tf:++c) 

{ 

8uax8q=8iimx8q+(ptil88  [a]  [c]  .shelll'^ptilaa  [a]  [c] .  ahelll) ; 
8UDy3q=8umysq-«- 

(atileaCa-t-l]  Cc]  .shelliaatileaCa+l]  [c]  .ahelll) ; 

> 

nauazaqsaumxaqeCV) ; 
naunyaq^aumyaqaCW) ; 

aiu8qz=aii]axe8UBZ : 
atuuqysaamyeaumy ; 

Bnmeratorsn8iuazy-(8iuuc*8amy) ; 

denoBlnators  (float)  aqrt  (  (double)  (naumzaq-aiuiiaqz)  e  (naunyaq-aunaqy)  ) ; 

if  (denoninator  o); 

{ 

pdel8Uffl«0 ; 

> 


/**  pdelaun  ia  the  correlation  coeff. 
pccnatch  ia  the  deaired  uinimun 
correlation 

tileratio  ia  the  rma  error  or 
difference  betveen  tilea 
tenpmg  ia  the  naximua  deviation 
between  tilea  or  acroaa  featureee/ 

if (denominator  !-  0); 

i 

pdelauBsnunerator/denoninator; 

> 

if  (pdel8un>>pccBatch  tileratio  <=  tenprng) 


/***  Tile  center  meeta  corr.  coeff,  map  full  tile  ***/ 
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K  the  reduced  set  sample  meets  the  conditional  statement  on  the  last  line,  the 
sample  is  remapped  in  a  similar  fashion  to  the  pattern  tile  routine. 

B.2.2  Search  tile  -  pattern  tile  matching 

The  correlation  for  the  expanded  set  is  computed  and  the  points  are  ranked  using 
a  routines  similar  to  the  ones  above.  The  following  criteria  are  applied  for  selection: 

if  (  pdel8um>-pccmatch  At  til«ratio<=tempmg 

AA  8lielllC0]>  8ma-t«Bipmg/2  AA  8helll[VW]<8m8-<-t«mprng/2) 

This  constrains  selection  to  a  specified  correlation,  maximum  change  in  temperature 
due  to  diurnal  heating,  and  a  maximum  temperature  range  for  the  tile. 

B.3  Vector  computations 

Computation  of  vector  magnitude  and  direction  is  a  three  step  process.  The  values 
for  the  search  tile  x  and  y  locations  are  subtracted  from  the  pattern  tile  values  and  length 
(magnitude)  is  computed.  Based  on  the  values  for  the  x  and  y  sums,  the  routine  selects 
a  specific  switch  subroutine  using  the  sign  of  the  sums  and  the  greater  of  their  absolute 
values. 

Vector  length  is  determined  by  a  simple  Pythagorean  routine: 

vconp  =  i-i2;  /*  x  and  y  movement  ♦/ 

hcomp  =  j-j2; 

length  =  (float) 8qrt(  (double)  (vcoB^evcomp  hcompehcomp) ) ; 

Vector  direction  (anglea)  is  computed  by  determining  the  angle  described  by  arccos  f. 
The  arrow  head  lines  issue  from  the  vector  end  at  45°  angles  from  the  line  (angleb,  anglec). 

anglea  -  (float)aco8 ( (double) (h  /  length)); 
anglebs  (pi/4)->'anglea ; 
anglecs (pi/4) -anglea ; 
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Direction  is  converted  to  degrees,  flow  rate  is  calculated,  azimuth  and  arrow  head 
endpoints  are  computed  based  on  the  region  in  which  the  vector  is  located  (Figure  A.2), 
and  the  image  graphic  is  initiated  with  a  function  call  as  shown  below: 

/*  conversion  for  radians  to  degrees  */ 

#deline  degconv  57.29577961 

vector [a] .anglea=anglaaedegconv: 
vectorCa] .flos«(length*pizaize)/dtg[0] .sec; 
vector  [a]  .angles  ((float)  (region  value)); 

/*♦  drawing  vector  **/ 
i  =  ptilesCa] [tilecenter] .i; 
il=  select [a]. i; 
j  =  ptiles [a] [tilecenter] . j ; 
jl=  select [a]. j; 

draB_line(image,zdim.ydiu, j ,i. jl.il, backvalue) ; 

/**  drawing  arrow  heads  **/ 

j2=j+icos; 

i2=i+isin; 

draw_line ( inage , zdin , ydin, j , i , j  2 , i2 . backvalue) ; 

j2=j+icoal; 

i2ai-isinl ; 

draw.line ( inage , zdin , ydin , j . i , j  2 , i2 , backvalue ) ; 


Appendix  C 


THERMAL  SIGNAL  MAPS 


This  appendix  contains  image  signal  three  dimensional  representations  generated 
using  MATLAB  for  the  mouth  of  the  Delaware  Bay.  These  figures  are  included  to  show 
the  signal  variations  that  are  used  by  the  SEDA  algorithm  to  select  features. 
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Cape  May,  New  Jersey 


Thermal  map:  arbitrary  scale.  Azimuth:  315 

Figure  C.2;  Thermal  signal  map,  Delaware  Bay  and  coastal  region,  0100  hours  local, 
11  march  1989.  Land  and  water  signals  are  nearly  Indistingmshable  due 
to  low  feature  temperatures  (5‘’  C),  which  are  approaching  the  lower  signal 
threshold  for  the  image. 
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Cape  May,  New  Jersey 


Thermal  map:  arbitrary  scale.  Azimuth:  315 

Figure  C.S:  Thermal  signal  map,  Delaware  Bay  and  coastal  region,  1158  hours  local, 
28  May  1989.  Smooth  ocean  signal  surface  indicates  formation  of  a  surface 
mixed  layer.  Elevated  land  signals  are  due  to  higher  incident  solar  flux  and 
seasonal  warming.  Extremely  low  values  on  left  edge  were  caused  by  clouds 
over  Lewes,  Delaware. 
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Cape  May,  New  Jersey 


Thermal  map:  arbitrary  scale.  Azimuth:  315 


Figure  C.4:  Thermal  signal  map,  Delaware  Bay  and  coastal  region,  0201  hours  local,  29 
May  1989.  Land  temperatures  lower  due  to  diurnal  heat  loss. 
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Cape  May,  New  Jersey 


Thermal  map:  arbitrary  scale,  Azimuth:  3 IS 

Figure  C.5:  Thermal  signal  map,  Delaware  Bay  and  coastal  region,  1303  hours  local,  10 
June  1989.  High  land  and  low  water  signals  due  to  diurnal  heating  from 
increased  incident  solar  flux. 
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Cape  May,  New  Jersey 


Thermal  map:  arbitrary  scale,  Azimuth:  315 

Figure  C.6:  Thermal  signal  map,  Delaware  Bay  and  coastal  region,  0117  hours  local,  11 
June  1989.  Land  values  are  lower  water  due  to  diurnal  heat  loss.  Delaware 
Bay  temperatures  (rear)  are  2-3°  C  higher  than  coastal  waters  (foreground). 
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Cape  May,  New  Jersey 


Thermal  map:  arbitrary  scale,  Azimuth:  315 

Figure  C.7:  Thermal  signal  map,  Delaware  Bay  and  coastal  region,  1308  hours  local,  29 
June  1989.  Signal  drop  in  foreground  due  to  clouds. 


Appendix  D 


TABULAR  DATA 


This  appendix  contains  tabular  output  data  for  one  program  run  in  the  period 
specified  in  the  Delaware  Bay  area.  Each  data  set  is  comprised  of  three  parts: 

•  Vector  direction  and  magnitude  information 

•  Search  tile  -  pattern  tile  correlation  information 

•  Program  run  parameters 


D.l  10-11  March  1080 


Froa:  89069170021.  To:  89069184165: 

Elapsed  time:  0  DayCs),  Total:  1  Hours,  41  Minutes,  34 
Search  box  diB:15  Offset:  10  Tile  dia;  11  X  11  luaber  of 
pt  Ingth(ka)  az(deg)  flow  a/s 


0 

0.00 

0.00 

0.00 

1 

0.00 

0.00 

0.00 

2 

0.00 

0.00 

0.00 

3 

0.00 

0.00 

0.00 

4 

0.00 

0.00 

0.00 

6 

0.00 

0.00 

0.00 

6 

0.00 

0.00 

0.00 

7 

0.00 

0.00 

0.00 

8 

0.00 

0.00 

0.00 

9 

0.00 

0.00 

0.00 

10 

0.00 

0.00 

0.00 

11 

0.00 

0.00 

0.00 

12 

0.00 

0.00 

0.00 

13 

0.00 

0.00 

0.00 

14 

0.00 

0.00 

0.00 

16 

0.00 

0.00 

0.00 

16 

0.00 

0.00 

0.00 

17 

0.00 

0.00 

0.00 

18 

0.00 

0.00 

0.00 

19 

0.00 

0.00 

0.00 

Seconds 
saaples  31 
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70 


20 

0.00 

0.00 

21 

0.00 

0.00 

22 

0.00 

0.00 

23 

0.00 

0.00 

24 

0.00 

0.00 

2S 

0.00 

0.00 

26 

0.00 

0.00 

27 

0.00 

0.00 

28 

0.00 

0.00 

29 

0.00 

0.00 

30 

0.00 

0.00 

Average  flov  distance: 


0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00  ka.  Avg  Az:  0.00  deg, 


Avg  IloB  rate;  0.00  m/s 


Search  tile  data  (z,7  coord  lor  center  subtile) 


ptile 

stile 

1  X 

y 

CC  val 

Pt  temp  St 

temp 

del  T 

Pt  sum 

St  sum 

Ratio 

0 

0 

137 

221 

0.46 

6.06 

4.26 

0.79 

4892 

4125 

1.19 

1 

0 

0 

0 

0.00 

0.00 

0.00 

0.00 

4793 

0 

0.00 

2 

0 

0 

0 

0.00 

0.00 

0.00 

0.00 

5249 

0 

0.00 

3 

0 

0 

0 

0.00 

0.00 

0.00 

0.00 

6318 

0 

0.00 

4 

0 

0 

0 

0.00 

0.00 

0.00 

0.00 

4818 

0 

0.00 

6 

0 

119 

243 

0.10 

4.97 

3.64 

1.44 

4816 

3424 

1.41 

6 

0 

130 

243 

0.12 

5.38 

4.56 

0.82 

5210 

4414 

1.18 

7 

0 

70 

264 

0.50 

6.50 

2.82 

2.68 

6327 

2732 

1.95 

8 

0 

109 

264 

0.30 

4.80 

3.48 

1.32 

4643 

3367 

1.38 

9 

0 

120 

254 

0.00 

4.89 

4.29 

0.60 

4730 

4152 

1.14 

10 

0 

79 

266 

0.32 

4.67 

3.04 

1.53 

4421 

2938 

1.60 

11 

0 

0 

0 

0.00 

0.00 

0.00 

0.00 

4604 

0 

0.00 

12 

0 

122 

266 

0.26 

4.87 

4.46 

0.42 

4718 

4316 

1.09 

13 

0 

0 

0 

0.00 

0.00 

0.00 

0.00 

6623 

0 

0.00 

14 

0 

87 

276 

0.20 

4.63 

3.47 

1.16 

4482 

3360 

1.33 

IS 

0 

0 

0 

0.00 

0.00 

0.00 

0.00 

4282 

0 

0.00 

16 

0 

116 

276 

0.16 

4.72 

4.06 

0.66 

4672 

3930 

1.16 

17 

0 

0 

0 

0.00 

0.00 

0.00 

0.00 

5086 

0 

0.00 

18 

0 

0 

0 

0.00 

0.00 

0.00 

0.00 

6526 

0 

0.00 

19 

0 

0 

0 

0.00 

0.00 

0.00 

0.00 

4349 

0 

0.00 

20 

0 

108 

287 

0.03 

4.60 

3.66 

0.94 

4366 

3460 

1.26 

21 

0 

120 

287 

0.02 

4.56 

4.01 

0.66 

4418 

3881 

1.14 

22 

0 

0 

0 

0.00 

0.00 

0.00 

0.00 

6379 

0 

0.00 

23 

0 

0 

0 

0.00 

0.00 

0.00 

0.00 

4176 

0 

0.00 

24 

0 

0 

0 

0.00 

0.00 

0.00 

0.00 

4273 

0 

0.00 

25 

0 

0 

0 

0.00 

0.00 

0.00 

0.00 

4876 

0 

0.00 

26 

0 

0 

0 

0.00 

0.00 

0.00 

0.00 

6402 

0 

0.00 

27 

0 

91 

309 

0.12 

4.31 

3.56 

0.74 

4168 

3447 

1.21 

28 

0 

104 

309 

0.25 

4.44 

3.76 

0.69 

4301 

3636 

1.18 

29 

0 

116 

309 

0.41 

4.91 

3.31 

1.59 

4761 

3208 

1.48 

30 

0 

0 

0 

0.00 

0.00 

0.00 

0.00 

6671 

0 

0.00 

System  sorking  values; 
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Inage  Xdlm:  612  Ydin:  612,  Pixsize:  1100 

Search  Area  xmin:  40  ymin:  200  xnaz:  140  ymax:  320,  Search  Option:  edge 

Ptile:  Edge  tile:  6,  Ptile:  11,  Offset:  10,  Grid  dim:  40,  Thresh:  6, 
Range:  12,  Naskl:  26,  I(ask2:  61 
Areal  SST:  6.62,  Area2  SST:  6.26 

Stile:  search  dim:  4,  offset:  1,  Allosed  RMS  temp  change:  4.00 
Alloved  temp  range  from  tile  mean:  2.000,  CC  match(desired) :  0.600, 

Flos  rate(desired) :  0.600,  Max  allosed  flos:  0.800 
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D.2  28-29  May  1989 


From:  89148166803.  To:  89149070111: 


Elapsed  time:  1  DayCs) ,  Total:  9  Hours.  67 
Search  box  dim: 103  Offset:  10  Tile  dim:  11  X 
pt  Ingth(km)  az(deg)  flow  m/s 


0 

48.27 

39.96 

0.92 

1 

69.46 

137.73 

1.14 

2 

61.74 

131.08 

0.99 

3 

0.00 

0.00 

0.00 

4 

0.00 

0.00 

0.00 

6 

23.00 

180.00 

0.44 

6 

0.00 

0.00 

0.00 

7 

49.04 

129.21 

0.94 

8 

67.49 

130.06 

1.10 

9 

0.00 

0.00 

0.00 

10 

0.00 

0.00 

0.00 

11 

41.19 

6.67 

0.79 

12 

62.80 

142.70 

1.01 

13 

0.00 

0.00 

0.00 

14 

0.00 

0.00 

0.00 

16 

0.00 

0.00 

0.00 

16 

0.00 

0.00 

0.00 

17 

16.62 

129.81 

0.30 

18 

0.00 

0.00 

0.00 

19 

0.00 

0.00 

0.00 

20 

47.80 

164.22 

0.92 

21 

0.00 

0.00 

0.00 

Minutes ,  8  Seconds 

11  Humber  of  samples  22 


Average  flow  distance:  13.73  km,  Avg  Az:  126.76  deg.  Avg  flow  rate;  0.26  m/s 


Search  tile  data  (x.y  coord  for  center  subtile) 


ptile  stile  z 

y 

CC  val 

Pt  temp  St  temp  del  T 

Pt  sum 

St  sum 

Ratio 

0 

81 

163 

196 

0.88 

12.01 

12.30  -0.30 

11621 

11909 

0.98 

1 

116 

103 

287 

0.72 

16.44 

14.68  1.86 

16916 

14110 

1.13 

2 

33 

106 

288 

0.63 

16.34 

14.69  0.66 

14847 

14222 

1.04 

3 

0 

119 

264 

0.02 

12.86 

13.71  -0.86 

12447 

13268 

0.94 

4 

0 

66 

266 

0.61 

14.96 

11.21  3.76 

14486 

10866 

1.33 

6 

14 

106 

288 

0.64 

12.31 

14.70  -2.39 

11916 

14231 

0.84 

6 

1 

126 

266 

0.28 

12.89 

13.17  -0.27 

12481 

12747 

0.98 

7 

44 

113 

307 

0.63 

13.40 

14.93  -1.63 

12972 

14460 

0.90 

8 

609 

160 

313 

0.70 

11.74 

13.69  -1.86 

11363 

13159 

0.86 

9 

0 

119 

276 

0.17 

12.30 

12.40  -0.09 

11908 

11999 

0.99 

10 

0 

136 

276 

0.24 

12.96 

13.42  -0.46 

12660 

12995 

0.97 

11 

19 

79 

246 

0.78 

14.19 

17.20  -3.01 

13740 

16661 

0.83 

12 

62 

122 

329 

0.64 

12.33 

13.96  -1.62 

11938 

13511 

0.88 

13 

0 

101 

287 

0.04 

12.46 

16.08  -3.63 

12066 

15566 

0.77 

14 

0 

117 

287 

0.16 

12.43 

12.60  -0.07 

12028 

12099 

0.99 

16 

0 

107 

298 

0.23 

13.03 

14.29  -1.26 

12614 

13829 

0.91 

16 

0 

126 

298 

0.17 

12.20 

13.18  -0.98 

11813 

12762 

0.93 

17 

63 

148 

308 

0.60 

12.30 

14.01  -1.72 

11902 

13663 

0.88 
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18 

0 

86 

300 

0.28 

13.80 

11.69 

2.20 

13444 

11316 

1.19 

19 

0 

98 

309 

0.04 

13.42 

9.98 

3.44 

12992 

9663 

1.34 

20 

716 

122 

356 

0.79 

13.16 

14.35 

-1.19 

12741 

13889 

0.92 

21 

0 

128 

309 

0.06 

11.87 

11.74 

0.13 

11490 

11366 

1.01 

Systea  vorking  values: 

Inage  Xdim:  512  Ydiin:  612.  Pixsize:  1100 

Search  Area  rain:  40  ymin:  200  zmax:  140  ymax:  320,  Search  Option:  edge 

Ptile:  Edge  tile:  6,  Ptile:  11,  Offset:  10,  Grid  dim;  40,  Thresh:  6, 
Range:  12,  Maskl:  86.  Mask2:  140 
Areal  SST:  16.60,  Area2  SST:  12.63 

Stile:  search  dim:  02,  offset:  1,  Alloved  RMS  temp  change:  4.00 
Allowed  temp  range  from  tile  mean:  2.000,  CC  match(desired) :  0.600, 

Flow  rate(desired) :  0.900,  Max  allowed  flow:  1.600 
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D.3  10-12  June  1089 

From;  89161180326.  To:  89162062806: 

Elapsed  time:  1  DayCs),  Total:  12  Hours,  26  Minutes,  19  Seconds 
Searcli  box  dim:83  Ollset:  10  Tile  dim:  11  X  11  lumber  of  s2UBples  19 
pt  Ingth(km)  az(deg)  flow  m/s 


0 

0.00 

0.00 

0.00 

1 

0.00 

0.00 

0.00 

2 

4.00 

180.00 

0.10 

3 

11.06 

96.19 

0.27 

4 

0.00 

0.00 

0.00 

6 

3.00 

180.00 

0.07 

6 

7.00 

180.00 

0.17 

7 

21.26 

48.81 

0.62 

8 

0.00 

0.00 

0.00 

9 

7.07 

171.87 

0.17 

10 

39.86 

197.63 

0.98 

11 

0.00 

0.00 

0.00 

12 

37.34 

172.30 

0.92 

13 

9.22 

319.40 

0.23 

14 

27.46 

169.61 

0.68 

16 

12.00 

180.00 

0.30 

16 

31.30 

26.67 

0.77 

17 

6.39 

201.80 

0.13 

18 

36.26 

24.44 

0.89 

Average  Iloe  distance:  3.77  km.  Avg  Az:  141.46  deg,  Avg  tlov  rate:  0.09  m/s 


Search  tile  data  (z.y  coord  lor  center  subtile) 


ptile  stile  z 

y 

CC  val 

Pt  temp  St  temp 

del  T 

Pt  sum  St  sum 

Ratio 

0 

0 

0 

0 

0.00 

0.00 

0.00 

0.00 

26302 

0 

0.00 

1 

0 

0 

0 

0.00 

0.00 

0.00 

0.00 

26832 

0 

0.00 

2 

163 

123 

236 

0.62 

18.70 

16.11 

2.68 

18099 

16699 

1.16 

3 

98 

147 

233 

0.81 

18.60 

16.34 

2.16 

17912 

16820 

1.13 

4 

0 

0 

0 

0.00 

0.00 

0.00 

0.00 

21696 

0 

0.00 

6 

136 

74 

246 

0.66 

22.12 

19.76 

2.38 

21414 

19116 

1.12 

6 

192 

116 

260 

0.64 

18.46 

16.44 

2.02 

17866 

16911 

1.12 

7 

8 

149 

229 

0.76 

19.03 

16.09 

2.94 

18423 

16676 

1.18 

8 

0 

0 

0 

0.00 

0.00 

0.00 

0.00 

23663 

0 

0.00 

9 

46 

66 

261 

0.66 

21.06 

18.36 

2.69 

20372 

17767 

1.16 

10 

669 

102 

292 

0.87 

17.80 

16.19 

2.61 

17231 

14701 

1.17 

11 

0 

0 

0 

0.00 

0.00 

0.00 

0.00 

23676 

0 

0.00 

12 

138 

114 

302 

0.69 

17.67 

17.60 

0.07 

17006 

16938 

1.00 

13 

116 

114 

268 

0.83 

18.67 

16.89 

2.67 

17971 

16383 

1.17 

14 

191 

111 

303 

0.74 

17.96 

17.20 

0.76 

17381 

16664 

1.04 

16 

199 

106 

299 

0.90 

18.38 

16.12 

2.26 

17796 

16607 

1.14 

16 

78 

107 

270 

0.66 

17.66 

16.63 

1.92 

16990 

16133 

1.12 

17 

161 

103 

303 

0.88 

18.67 

16.73 

2.84 

17973 

16226 

1.18 

18 

16 

107 

276 

0.69 

17.92 

16.24 

1.68 

17361 

16723 

1.10 

System  vorking  values: 
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laag*  Zdia:  612  Ydim:  612,  Pixsiza:  1100 

Search  Area  zmin:  40  200  zmax:  140  ymax:  320,  Search  Option:  edge 

Ptile:  Edge  tile:  6,  Ptile:  11,  Offset:  10,  Grid  dim:  40,  Thresh:  6, 
Range:  12.  Naskl:  132,  Nask2:  182 
Areal  SST:  21.26,  Area2  SST:  18.00 

Stile:  search  dim:  72,  offset:  1,  Aliened  RMS  temp  change:  3.00 
Aliened  temp  range  from  tile  mean:  1.600,  CC  matchCdesired) :  0.600, 

Flon  rata(dasired) :  0.900,  Max  aliened  flon:  1.600 
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D.4  20-30  June  1080 


Froa:  89180180847.  To:  89181063266: 


Elapsod  tiaa:  1  Day (a).  Total:  12  Hours,  24  Minutes,  9  Seconds 
Search  box  din: 83  OTlset:  10  Tile  din:  11  X  11  lunber  of  samples  26 
pt  Ingth(km)  az(deg)  flow  n/s 


0 

0.00 

0.00 

1 

0.00 

0.00 

2 

0.00 

0.00 

3 

23.71 

162.36 

4 

23.41 

160.02 

6 

43.14 

46.94 

6 

0.00 

0.00 

7 

41.11 

48.96 

8 

0.00 

0.00 

9 

1.00 

180.00 

10 

0.00 

0.00 

11 

46.80 

63.88 

12 

0.00 

0.00 

13 

24.76 

46.64 

14 

2.24 

163.43 

16 

0.00 

0.00 

16 

0.00 

0.00 

17 

60.46 

60.63 

18 

19.31 

168.76 

19 

0.00 

0.00 

20 

4.00 

360.00 

21 

6.00 

216.87 

22 

21.00 

360.00 

23 

23.77 

337.76 

24 

10.00 

143.13 

26 

26.02 

367.80 

Average  flow 

distance 

0.00 

0.00 

0.00 

0.68 

0.68 

1.06 

0.00 

1.01 

0.00 

0.02 

0.00 

1.13 

0.00 

0.61 

0.06 

0.00 

0.00 

1.24 

0.48 

0.00 

0.10 

0.12 

0.62 

0.69 

0.26 

0.64 

8.10  km,  Avg  Az:  63.66  deg. 


Avg  lloB  rate:  0.20  m/s 


Search  tile  data  (x,y  coord  for  center  subtile) 


ptile  stile  z 

y 

CC  val 

Pt  temp  St  temp 

del  T 

Pt  sum 

St  sum 

Ratio 

0 

14 

66 

210 

0.24 

27.40 

22.67 

4.73 

26628 

21949 

1.21 

1 

99 

81 

210 

0.27 

28.61 

24.08 

4.63 

27695 

23312 

1.19 

2 

109 

119 

210 

0.06 

26.12 

21.22 

4.90 

26284 

20644 

1.23 

3 

391 

72 

253 

0.82 

24.42 

21.64 

2.88 

23642 

20861 

1.13 

4 

683 

76 

266 

0.82 

22.87 

20.02 

2.85 

22136 

19376 

1.14 

6 

35 

137 

213 

0.63 

22.02 

18.90 

3.12 

21319 

18295 

1.17 

6 

0 

118 

243 

0.06 

21.36 

18.83 

2.64 

20679 

18224 

1.13 

7 

63 

168 

216 

0.68 

20.81 

21.00 

-0.19 

20143 

20327 

0.99 

8 

0 

60 

264 

0.41 

26.61 

22.64 

3.06 

24789 

21823 

1.14 

9 

484 

73 

266 

0.90 

21.83 

21.43 

0.40 

21133 

20742 

1.02 

10 

22 

92 

264 

0.17 

27.02 

22.22 

4.80 

26151 

21509 

1.22 

11 

79 

145 

227 

0.64 

20.61 

20.29 

0.22 

19865 

19641 

1.01 

12 

1 

120 

264 

0.11 

21.14 

18.63 

2.61 

20468 

17938 

1.14 

13 

101 

76 

248 

0.86 

24.69 

22.86 

1.74 

23806 

22124 

1.08 

77 


14 

426 

76 

267 

0.89 

20.40 

19.60 

0.80 

19762 

18976 

1.04 

16 

2 

91 

266 

0.19 

20.34 

16.33 

4.01 

19687 

15808 

1.26 

16 

0 

104 

266 

0.06 

20.17 

19.97 

0.20 

19629 

19332 

1.01 

17 

3 

163 

233 

0.67 

21.04 

21.90 

-0.86 

20371 

21197 

0.96 

18 

396 

88 

294 

0.84 

18.27 

19.23 

-0.96 

17681 

18614 

0.96 

19 

0 

70 

287 

0.06 

26.82 

22.80 

4.03 

26966 

22067 

1.18 

20 

4 

82 

283 

0.80 

18.20 

18.16 

0.06 

17621 

17666 

1.00 

21 

129 

91 

291 

0.92 

20.24 

19.90 

0.34 

19692 

19262 

1.02 

22 

160 

82 

277 

0.70 

18.13 

18.11 

0.01 

17646 

17534 

1.00 

23 

100 

129 

276 

0.61 

21.83 

21.53 

0.31 

21136 

20837 

1.01 

24 

13 

87 

317 

0.63 

18.63 

19.48 

-0.96 

17934 

18869 

0.95 

26 

4 

91 

283 

0.84 

20.56 

19.63 

0.93 

19904 

18999 

1.05 

System  vorking  values: 

Image  Xdim:  612  Ydim:  612,  Pixsize:  1100 

Seuch  Area  xmin:  40  ymin:  200  zmaz:  140  ymaz:  320,  Search  Option:  edge 

Ptlle:  Edge  tile:  6,  Ptile:  11,  Offset:  10,  Grid  dim:  40,  Thresh:  6, 
Range:  12,  Naskl:  124,  Mask2:  196 
Areal  SST:  22.00,  Area2  SST:  18.00 

Stile:  search  dim:  72,  offset:  1,  Alloved  RMS  temp  change;  6.00 
Aliened  temp  range  from  tile  mean:  2.500,  CC  match(desired) :  0.600, 

Flow  rate(dasired) :  0.900,  Max  aliened  flon:  1.600 
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